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SUMMARY  OF  SCIENTIFIC  PROGRESS  TO  DATE 


Here  brief  summaries  are  given  for  topics  where  significant  progress 
has  been  achieved,  but  where  as  yet  no  papers  have  been  prepared  for 
publication, 

1.  Non-Linear  Theory  (I.  B.  Bernstein  and  L.  Friedland) 

The  non-linear  theory  of  a  steady  state  free  electron  laser  amplifier 
consisting  of  a  relativistic  electron  beam  in  helical  pump  and  axial  con¬ 
fining  magnetic  fields  illuminated  by  monochromatic  radiation  has  been  form¬ 
ulated,  The  introduction  of  Lagrangian  variables  permits  an  exact  integra¬ 
tion  of  the  continuity  equation  to  express  the  current  density  at  a  general 
field  point  in  terms  of  the  current  density  at  injection  and  the  velocity  at 
the  field  point  in  question*  The  axial  electric  field  is  found  in  a  similar 
way  by  integrating  Poisson’s  equation.  The  parallel  momentum  equation  is 
reduced  to  an  ordinary  differential  equation  and  permits  the  electron  orbits 
to  turn  in  the  wave  frame,  corresponding  as  one  advances  along  the  beam  in 
the  direction  of  propagation  to  the  development  of  multistreaming*  The 
reduced  equations,  and  the  coupled  wave  equations  for  the  transverse  vector 
potential  are  being  analyzed,  and  suitable  numerical  methods  for  their  solu¬ 
tion  are  being  determined.  Preliminary  numerical  studies  have  begun. 

2.  Quasi-L inear  Theory  (I.  B.  Bernstein  and  L.  Friedland) 

A  quasi-linear  theory  of  the  steady  free  electron  laser  amplifier  con¬ 
sisting  of  a  helical  pump  and  relativistic  electron  beam  illuminated  by  a 
monochromatic  wave  has  been  formulated.  A  linearized  WKB  description  is 
employed  to  describe  the  electromagnetic  field,  and  the  steady  state  density, 
mean  velocity,  etc.  are  determined  correct  to  second  order  in  the  high  fre¬ 
quency  amplitude.  The  helical  pitch  is  allowed  to  vary  with  position  so  as 
to  keep  the  system  in  a  configuration  of  maximum  gain  as  the  beam  deceler¬ 
ates  on  converting  its  energy  into  radiation.  The  analysis  of  the  results 
is  continuing,  numerical  studies  are  underway,  and  a  preliminary  manuscript 
has  been  prepared. 

3.  Exact  Magnetic  Fields  of  a  Bifilar  Helix  (S.  Y.  Park,  R.  A,  Smith,  and 
J.  L.  Hirshf ield) 

The  static  magnetic  field  of  a  bifilar  helix  has  been  derived  exactly. 
The  components  B]_,  B2,  and  B3  have  each  been  determined  as  functions  of  r, 
and  0-kz,  where  2nYk  is  the  period  of  the  double  helix,  and  where  (e^,  §2* 

£3)  is  the  basis  vector  set  dictated  by  helix  symmetry.  In  usual  FEL  theory 
B-j^  -  B3  *  0,  and  B2  *  const.  The  theory  developed  here  gives,  for  a  helix 
of  specified  radius  and  period,  the  radial  beam  positions  for  which  the 
usual  assumption  is  valid.  Higher  order  multi-pole  fields  are  significant 
off  axis,  where  annular  beams  are  commonly  placed. 

4.  Non-Helical  Orbits  in  a  Free  Electron  Laser  (S.  Y.  Park,  R.  A.  Smith  and 
J.  L.  Hirshfield) 

An  exact  analytic  solution  has  been  found  for  the  orbit  of  a  charged 
particle  moving  in  the  combined  magnetic  field  of  a  bifilar  helical  wiggler 
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and  a  uniform  solenoid.  The  axial  velocity  is  shown  to  satisfy  an  anharmonic 
oscillator  equation  whose  solutions  are  Jacobian  elliptic  functions.  Other 
velocity  components  are  readily  determined  from  the  axial  velocity.  Solu¬ 
tions  are  classified  according  to  the  type  of  Jacobian  elliptic  function, 
which  in  turn  is  determined  by  magnetic  field  parameters  and  initial  particle 
conditions.  This  exact  solution  is  to  be  used  to  calculate  the  spectrum  of 
single-particle  spontaneous  emission;  significant  emission  at  harmonic  wave¬ 
lengths  is  indicated. 

5.  Electron  Beam  Analyzer  (P.  Avivi,  F.  Dothan,  A.  Fruchtman,  and  J.  L, 
Hirshfield) 

Electromagnetic  gain  in  a  FEL  is  a  sensitive  function  of  the  momentum 
distribution  of  the  beam  electrons.  The  physical  mechanism  which  dominates 
the  gain  process  may  even  change  as  one  goes  from  a  "cold"  beam  to  a  beam 
with  a  finite  skewed  spread  in  momentum.  It  is  thus  crucial  to  measure  the 
momentum  distribution  of  the  electron  beam  in  a  FEL  in  order  to  interpret 
the  results.  With  this  in  mind  we  have  designed  a  90°  magnetic  deflection 
analyzer  which  is  capable  of  determining  average  electron  momentum  and  rela¬ 
tive  distribution  of  momentum,  for  a  small  sample  of  the  beam.  At  a  deflec¬ 
tion  angle  of  90°  the  analyzer’s  dispersion  is  a  maximum;  momentum  resolu¬ 
tion  of  0.1%  appears  possible.  This  analyzer  is  now  under  construction,  in 
two  stages.  In  the  first  stage,  about  to  be  installed,  the  analyzer  responds 
in  a  single  channel  (i.e.,  in  a  single  momentum  bin)  on  each  accelerator 
shot.  In  the  second  stage,  a  multichannel  collector  and  readout  will  be 
possible  on  each  shot. 

6.  Operation  of  Febetron  Accelerator  (A.  Ljudmirsky,  J.  Smulovits,  F.  Dothan, 
H.  Sharon,  and  J.  L.  Hirshfield) 

The  ^  400  keV  Febetron  accelerator  has  been  successfully  operated  in 
the  magnetic  field  and  vacuum  tank  in  which  the  free  electron  laser  experi¬ 
ments  are  to  be  conducted.  The  manufacturer’s  Blumlein  pulse  line  has  been 
removed  so  that  the  50  nsec  Marks  bank  is  terminated  directly  by  the  foil¬ 
less  field  emission  diode.  Diode  impedance  is  adjusted,  both  empirically 
and  with  the  guidance  of  the  NRL  experience,  by  shaping  the  cathode  stem 
and  the  anode  geometry.  Stable  beams  in  the  current  range  100-1000  A  have 
been  successfully  generated  and  launched  along  the  full  axis  of  the  combined 
guide  and  wiggler  magnetic  fields  of  the  device.  Electromagnetic  radiation 
and  gain  measurements  at  35  GHz  are  about  to  begin. 
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Free-Electron  Laser  with  a  Strong  Axial  Magnetic  Field 

L.  Friedland  and  J.  L.  Hirshfield 
Department  of  Engineering  and  Applied  Science,  Mason  Laboratory , 

Yale  University ,  New  Haven,  Connecticut  06520 
(Received  27  February  1980) 

A  small-signal  theory  is  given  for  gain  in  a  free-electron  laser  comprising  a  cold 
relativistic  electron  beam  in  a  helical  periodic  transverse,  and  a  strong  uniform  axial, 
magnetic  Held.  Exact  finite- amplitude,  steady-state  helical  orbits  are  included.  If 
perturbed,  these  orbits  oscillate  about  equilibrium,  so  that  substantial  gain  enhancement 
can  occur  if  the  electromagnetic  perturbations  resonate  with  these  oscillations.  This 
gain  enhancement  need  not  be  at  th^  cost  of  frequency  upshift. 

PACS  numbers:  42.55.-f,  41.70.+t,  41.80.Dd 

known  result®  for  G0  in  the  single-particle  limit 
(i.e.,  when  collective  effects  are  negligible) 

G0  * (u>,  C /V)2 (fco<L/2y)3F' '(9).  (1) 

Here  and  y  are  the  beam  plasma  frequency 
Ne2/m€0  and  normalized  energy  W /me2,  k0  and  l 
are  the  heitcal  transverse  magnetic  field  wave 
number  2 tt /l  and  normalized  strength  eBjmck^, 

L  is  the  interaction  length,  and /'(l)  «(ti/tf0)(sin0/ 
@)2  is  the  line-shape  factor,  with  9 
-v30/c)](L/2c)f  where  t/30  is  the  unperturbed  elec¬ 
tron  axial  velocity.  The  peak  gain  occurs  at  6 
*1.3,  where  F'(9)  =0.54.  For  example,  withy 


Intensive  activity  is  underway  to  exploit  the 
gain  properties  of  a  relativistic  electron  beam 
undulating  in  a  periodic  transverse  magnetic  field. 
Such  free-electron  laser  (FEL)  configurations 
have  provided  oscillation  at  3.4  (Ref.  1)  and  400 
jim,2  and  amplification  at  10.6  jim.3  Theory  has 
advanced  apace,4  and  elaborate  schemes  have 
been  proposed  for  obtaining  high  FEL  efficiency.5 
A  factor  which  limits  the  practical  application  of 
this  interaction  at  wavelengths  shorter  than  per¬ 
haps  a  few  microns  is  the  rapid  decrease  in 
small- signal  gain  G0  as  the  electron  energy  in¬ 
creases.  This  is  shown  explicitly  in  the  well- 
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=  10,  l  =  1.05  cm,  w,  =  5xlo7  sec’1,  4=1  (^  =  10.2 
kG),  and  L  =130  cm,  the  peak  gain  is  G^* 0.002 47 
at  a  wavelength  of  105  jim.  For  y  =100,  Z=10.5 
cm,  up  =2xl09  sec"1,  4=1  (Bx  =  l,02  kG),  and  L 
=  260  cm,  the  peak  gain  is  =0.003  16  at  a  wave¬ 
length  of  10.5  Mm,  These  gain  values  may  be 
large  enough  to  sustain  oscillations  if  highly  re¬ 
flecting  mirrors  are  judiciously  added  but  the 
strong  helical  fields  required  (particularly  the 
10.2-kG  case)  may  be  beyond  the  capability  of 
present  superconducting  coil  technology.7 

A  suggestion  has  appeared  for  enhancing  the 
small-signal  gain  above  values  given  by  Eq.  (1) 

(or  for  achieving  comparable  gains  with  smaller 
B  J  by  employing  a  strong  axial  magnetic  field  so 
as  to  exploit  resonance  between  the  cyclotron  fre¬ 
quency  and  the  undulatory  frequency.8  The  pres¬ 
ent  Letter  presents  a  single-particle  derivation 
for  the  small-signal  gain  of  a  FEL  in  a  uniform 
axial  magnetic  field  B0.  We  shall  demonstrate 
that  careful  adjustment  of  the  system  parameters 
will  allow  enhancement  of  the  FEL  small- signal 
gain  by  an  order  of  magnitude  or  more  (for  the 
above  examples)  without  increasing  the  undulatory 
velocity.  This  result  goes  beyond  that  predicted 
by  Sprangle  and  Granatstein®  who  have  suggested 
that  the  only  effect  of  the  axial  magnetic  field 
would  be  to  add  a  multiplicative  factor  (1  -Q/ 
k0 cy)’2  to  Eq.  (1),  due  to  the  aforementioned  res¬ 
onance  giving  an  enhanced  undulatory  velocity  vXf 
where  f 1  -eBjm.  This  result  is  in  fact  predicted 
by  our  analysis  as  a  limiting  case.  Of  course, 
any  mechanism  which  increases  the  undulatory 
velocity  v±  would  increase  the  gain,  but  this 
would  also  reduce  the  relativistic  frequency  up¬ 
shift,  since 

^^(l  -u30/c)‘l  =  2y2*0c(l  +y2i>x2/c2)’1 . 

If,  for  example,  yv  Jc  =1  without  the  axial  mag¬ 
netic  field,  then  a  given  gain  enhancement  rj 
achieved  through  this  resonance  alone  would  re¬ 
sult  in  a  reduction  in  frequency  upshift  by  a  fac¬ 
tor  (1  +tj)/2.  The  process  we  describe  in  this 
Letter  will  be  shown  to  permit  significant  gain 
enhancement  without  undue  sacrifice  in  frequency 
upshift.  The  gain  enhancement  originates  when 
the  electromagnetic  perturbations  resonate  with 
the  natural  frequency  of  oscillation  of  electrons 
on  finite  amplitude  equilibrium  helical  orbits. 

Prior  workers  have  not  considered  this  effect. 

A  full  derivation  of  our  result  will  be  presented 
elsewhere.9  Exact  unperturbed  relativistic  orbits 
are  considered  in  the  customary  FEL  model  mag¬ 


netic  field 

B(r)  =B0e g  +Bx(ez  cos*^  +ey  sinfe^).  (2) 

These  orbits,  which  have  been  the  subject  of  re¬ 
cent  study,10  can  possess  more  than  one  steady 
state,  depending  upon  y ,  B0,  Bx,  and*0.  These 
steady  states  are  characterized  by  the  normalized 
velocity  components  (i.e.,  ut  = vjc ) 

^10  *^2Q  "■  Sl/c)  f 

(3) 

ttS0  =  (l-tt302  -y*2)1'2, 

where  the  basis  vectors  eL(z)  =  -e,  sinfe^  +e, 
xcoskgZ,  e2U)  s-lyCos^^ -eysinfe^,  ande3fc) 
=es  have  been  introduced  to  track  the  symmetry 
of  the  transverse  magnetic  field.  Figure  1  shows 
u30  vs  Sl/c  for  *o=6.0  cm"1,  4=1.0,  andy=10. 

For  G>flCrs^oct()/2"1)1/3“42/3]3/2>  is  seen 
that  only  one  branch  exists  (branch  C).  But  for 
Q  <  Hcr  two  additional  branches  (A  and  B)  are  al¬ 
lowed:  Branch  B  has  been  shown  to  be  unstable, 
in  that  the  orbits  exhibit  nonhelical,  highly  an- 
harmonic  motions,  while  branches  A  and  C  have 
orderly  helical  orbits.  Stability  is  insured  if  m2 
=a2-bd> 0,  where  a  =k0cu30l/yu20,  b  =au20/yu20, 
andd  =k0cZ/y.  The  quantity  ^  is  the  natural  res- 
onance  frequency  in  response  to  small  perturba¬ 
tions  of  the  orbit:  We  shall  show  that  strong  res¬ 
onance  response  of  the  electrons  to  electromag¬ 
netic  perturbation  can  lead  to  enhanced  FEL  gain 
for  small  m,  i.e.,  for  Q  close  to  £lcr. 

The  derivation  of  FEL  gain  proceeds  by  solving 
the  single-particle  equations  of  motion,  subject 
to  weak  electromagnetic  perturbing  fields  E 
=erB0  cos(fe*  -u>f)  and  B  =ey(*c/w)£0cos(*£  -  ut), 


SI/ c  (cm-1) 

FIG.  1.  Steady- state  normalized  axial  velocity  u  «,0 
as  a  function  of  normalized  axial  magnetic  field  ft/c. 
For  this  example  =  6.0  cm' l,  4  ==  1.0,  and  y  =  10. 
Gain  enhancement  discussed  in  this  work  is  for  orbits 
on  either  branch  A  or  branch  C . 
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about  the  equilibrium  orbits  on  either  branch  A  or  C  as  discussed  above.  These  equations  are 

-  (feoCw 3  -n/y)u,  -  (fe0c{/y)«3  -  (y/y)Wi  +  (e£l/mcy)(fccu3/u>  -  1),  (4) 

ua«-  (fe0cu3)  -O/y)^  -  (y/y)ua  +  fe£a/mcy)(*cu3/w  -  1),  (5) 

it2*(J*<fii/y)ux  +  (fee/u>  -u3)(y/y),  (6) 

where  y  *-  (E/mc)(uI£1  +ua£a)  and 

2(Ea+t£l)  *-£0exp{i[(fe0+*)u3Ct  -wf  +a]} 

with  a  the  random  initial  electron  phase.  When  time  variations  and  electromagnetic  fields  are  absent, 
Eqs.  (4)-(6)  lead  to  the  exact  steady  states  given  by  Eq.  (3).  To  linearize  Eqs.  (4)-(6),  we  introduce 
the  velocity  perturbations  w{  =u,  -u,0«al0  and  retain  only  the  lowest-order  quantities.  This  results  in 
wx  ^^2wl  *A£0cos@t  +a),  or 

wx  T^Sa  [cosf&t  +a)  -cos^f  cosa  +  (8/m)  sinjif  sina]  +  m’  lwx(0)  sinMf,  (7) 


where 


A  *(a  +0)(1  -tt30)+huao,  ^  =c(fe+fe0)«30- iil(0)=fe£o/2>7nc)(l-U30)sina, 
and  ^(0)  -0.  The  other  components  follow  from 

w2  +  fe£0/2mcy)(l  -«30-u202)  cos^f  +  a),  w2(0)  =0;  (8) 

and 


+(eE0/2mcy)u20(l  -m30)  cos(0f  +  a),  «/3(0)  =0.  (9) 

Equation  (7)  for  wx  exhibits  the  aforementioned  natural  resonance  at  frequency  m.  while  the  electromag¬ 
netic  perturbation  drives  the  transverse  motion  at  frequency  *3.  Gain  enhancement  can  be  expected 
when  n  is  close  to0. 

The  energy  gain  for  an  electron  is  calculated  from  (;nc/e)dy/dt =*  -WiEl0-w2E2Q-u20E2i.  The  first- 
order  variation  in  electric  field  £ai  originates  from  small  phase  variations  as  u3  changes.  Thus  this 
becomes 

(jnc/e)dy/dt  -m/2£20-  i£0(fe  +ko)cu2osin(0t  +a)fQ'dt'w:i{r).  (10) 

The  third  term  in  Eq.  (10)  is  much  larger  than  the  other  two  on  account  of  the  factor  k  +fc0.  The  domi¬ 
nant  single-particle  energy  transfer  in  the  FEL  (even  with  an  axial  magnetic  field)  is  seen  to  be  by 
work  ecuj£2  done  along  the  transverse  undulatory  motion,  enhanced  by  the  strong  variation  in  £2  as  its 
phase  varies  through  k>3.  The  energy  variation  [Eq.  (10)]  is  averaged  over  random  phase  a  to  give 
< dy/dt ),  which  in  turn  leads  to  the  gain  through  G  =  2(e0£02)"  lNmc2j*dt(dy/dt),  where  N  is  the  beam 
electron  density  and  T  *L/c  is  the  total  interaction  time  for  the  electrons  in  a  system  of  length  L. 

The  final  result  is 


F{d  ^)-f(e-(^) 

2  <p 


F(9+<?)-F(fl-<?)  _a__  r„,/nV  P(0+<p)-P(0-<p)f! 

"m*tL p{6) - J  ’ 


2  cp 


(ID 


where  0=0T/2,  <p=nT/ 2,  £(*)  =  ( sinx/x)2,  and 
PM  =xFM/ 2;  and  where  we  have  approximated 
(fc  +fc0)(l  -u30)=*fc0.  We  shall  examine  Eq.  (11)  in 
several  limits. 

For  only  the  terms  involving  F'(S)  and 

P'(0)  in  Eq.  (11)  are  significant,  and  on  branch 
A  the  latter  of  these  is  smaller  than  the  former 
by  at  lea*!  a  factor  2<p.  Thus  to  a  good  approxi¬ 


mation  we  may  write 

G( m  »  3)  =  Z{vi. pz/l&y)k0cu202T3F'(9) ,  (12) 

where  Z-  2  +  +  -“3c))*1]*  In  the  case 

where  the  axial  magnetic  field  is  absent,  fl  =  0, 

M  =  a  =  k0cu20  »  0,  and  u20  =  £/y.  Thus,  Za2  and 
Eq.  (12)  goes  over  to  Eq.  (1).  When  fl*0  and  ^ 
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»3,  gain  enhancement  can  be  achieved  as  claimed 
by  the  prior  workers,9  due  to  resonant  enhance¬ 
ment  of  iUq,  but  not  without  sacrificing  frequency 
upshift,  as  discussed  above. 

However  a  more  attractive  possibility  exists 
when  ji  is  small,  and  approaches  3.  Here  one 
can  approximate  Z  =*  ^2bd(\  -  u30)~l  »  1;  this  re¬ 
sults  from  resonance  between  the  electromagnetic 
perturbation  which  gives  oscillatory  motion  to  the 
electron  at  a  frequency  3,  close  to  its  natural 
oscillation  frequency  p.  Gain  enhancement  due 
to  large  Z  is  seen  to  be  possible  without  simul¬ 
taneously  increasing  u20,  so  that  the  desirable 
frequency  upshift  property  of  the  FEL  need  not 
be  sacrificed. 

We  define  a  gain  enhancement  factor  77  =  G/G0 
to  compare  two  free-electron  lasers,  identical 
except  that  one  has  a  strong  axial  magnetic  field, 
while  the  second  does  not.  In  the  first  laser,  the 
transverse  magnetic  field  B  L  is  reduced  so  that 
u20  is  the  same  for  both  lasers.  (This  assures 
that  both  enjoy  the  same  frequency  upshift.)  Then 

v  =  z{  1  -  [/X e  +  <p)  -  F{ 9  -  <p)  ]/2( PF'{ e)} .  (13) 

We  have  evaluated  Eq.  (13)  for  two  examples  with 
the  parameters  cited  in  the  first  paragraph  of  this 
Letter,  holding  ( fc?|  =  1.3  where  |/v(0)|  has  its 
maximum  value.  The  results  are  shown  in  Fig.  2 
for  the  y-  10  example.  In  Fig.  2(a)  we  plot  the 
gain  enhancement  factor  77  as  a  function  of  the 
transverse  magnetic  field  normalized  strength  £ 
for  the  FEL  with  the  axial  guide  magnetic  field. 
The  solid  curves  are  for  steady-state  orbits  on 
branch  C;  the  dashed  curves  for  branch  A .  On 
branch  A,  gain  occurs  for  0>O,  while  on  branch 
C  gain  occurs  for  0<O.  Two  transverse  magnetic 
fields  for  the  FEL  without  axial  field  correspond¬ 
ing  to  £,  =  1  and  0.5  are  shown.  Figure  2(b)  shows 
the  required  values  of  axial  guide  field.  One  sees 
a  gain  enhancement  of  31  (on  branch  C)  at  4  =  5 
xlO"3  for  an  axial  guide  field  of  102  kG.  The 
transverse  magnetic  field  required  is  reduced  to 
51  G,  and  the  gain  is  increased  to  0,0766  at  \ 

-  105  4m.  Higher  gain  is  predicted  on  branch  A. 
For  the  y  -  100  example  at  \  =  10.5  pm,  we  find  a 
gain  enhancement  of  16  (on  branch  A)  at  £  =  3xl0"2 
for  an  axial  guide  field  of  99.5  kG.  The  trans¬ 
verse  magnetic  field  required  is  reduced  to  20.6 
G,  and  the  gain  is  increased  to  0.0506. 

Of  course  when  the  predicted  single-pass  gain 
is  large  (say  >0.1)  this  theory  must  be  modified. 
Furthermore,  finite  electron  momentum  spread 
(neglected  here)  will  mitigate  against  gain,  as 
for  a  FEL  without  a  guide  field.  These  effects 


FIG.  2.  (a)  Gain  enhancement  I  J?j  and  (b)  correspond¬ 
ing  normalized  axial  magnetic  field  ft/c  ,  vs  transverse 
magnetic  field  parameter  4.  The  values  4o  ~  0.5  and 
1.0  are  for  the  FEL  without  axial  field,  and  provide 
the  same  u10  as  do  the  indicated  (smaller)  values  of 
4  for  the  FEL  with  the  indicated  axial  field  strength. 
Example  is  for  y~  10,  k0  =  6.0  cm’ l,  and  L  -  130  cm. 
Solid  curves,  orbits  on  branch  C;  dashed  curves, 
orbits  on  branch  A.  For  high  enhancement  values, 
such  as  on  the  £0  -  1-0  branch  A  example,  the  numeri¬ 
cal  precision  required  to  compute  accurate  results 
suggests  that  the  phenomenon  is  very  sensitive  to  the 
system  parameters . 


deserve  careful  study.  However,  to  the  extent 
that  these  effects  are  negligible,  our  theory 
shows  that  provision  of  a  strong  uniform  axial 
magnetic  field  can  allow  significant  small-signal 
gain  enhancement,  and  significant  reduction  in  the 
required  transverse  magnetic  field  strength  in  a 
FEL,  without  undue  compromise  in  operating  fre¬ 
quency  below  that  given  by  the  idealized  upshift 
value  2 y2k0c. 
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The  propagation  of  a  cold  relativistic  electron  beam  in  a  free  electron  laser  with  an  axial  guide  magnetic  field 
is  considered  The  possibility  of  several  steady-state  helical  trajectories  for  the  electrons  is  shown,  and  the 
stability  against  perturbations  and  accessibility  of  such  steady  states  is  considered.  Necessary  and  sufficient 
conditions  for  the  stability  are  derived  and  indicate  the  importance  of  the  transition  region  at  the  entrance  of 
the  laser.  Possible  modes  of  operation  of  the  laser  in  different  steady-state  regimes  are  suggested  and 
illustrated  by  numerical  examples. 


I.  INTRODUCTION 

The  propagation  of  a  relativistic  electron  beam  in 
transverse  periodic  magnetic  structures  has  been  stu¬ 
died  extensively  in  recent  years.  These  studies  were 
stimulated  dv  the  first  experimentally  successful  free 
electron  laser1  which  confirmed  the  theoretically  pre¬ 
dicted2  possibility  of  using  the  energy  stored  in  the 
relativistic  beam  as  a  source  of  short  wavelength  co¬ 
herent  radiation. 

The  most  frequently  used  periodic  magnetic  pump  field 
in  a  free  electron  laser  is  the  transverse  field  produced 
on  the  axis  of  a  double  helical  current  winding  with  equal 
and  opposite  currents  in  each  helix  (a  device  usually 
referred  to  as  a  magnetic  wiggler).  The  unperturbed 
motion  of  the  electrons  of  the  beam  in  the  wiggler  is 
quite  simple.  The  reason  for  simplicity  is  that  the 
magnetic  field  on  the  axis  of  a  wiggler  can  be  approxi¬ 
mately  described  by  a  transverse  vector  potential  A^z), 
depending  only  on  the  distance  z  along  the  axis.  There¬ 
fore,  the  canonical  transverse  momentum  Pi  =  y»ivA 

-  (e/c) Aa  of  an  electron  is  a  constant  of  motion,3  which 
with  the  conservation  of  energy  y  =  [l  -  ( vjc )2 

-  (t’n/c)2]"*72  =  const,  uniquely  defines  the  perpendicular 
and  parallel  components  vA  and  v„  of  the  velocity  of  the 
electrons  in  the  beam,  for  a  given  assignment  of  A^z)- 
It  can  now  be  easily  shown4  that  the  electrons  in  a  mag¬ 
netic  wiggler  have  helical  trajectories  with  the  same 
period  as  that  of  the  wiggler.  This  simple  model  of  the 
motion  has  been  exploited  in  many  theoretical  studies, 
describing  the  operation  and  parametric  behavior  of  the 
free  electron  laser.3 

In  all  experiments,  however,  there  is  also  an  axial 
guide  magnetic  field.1*6*7  This,  of  course,  increases 
the  number  of  parameters  characterizing  the  free  elec¬ 
tron  laser,  but  at  the  same  time  introduces  greater 
complexity  into  the  theory.  The  vector  potential  is  now 
dependent  on  x  and  y  and  the  perpendicular  canonical 
momentum  Px  is  no  longer  a  constant  of  motion;  as  a 
result,  in  general,  no  simple  analytic  solution  for  the 
electron  trajectories  can  be  found.  Although  so  called 
“steady- state"  helical  trajectories  with  the  same  period 
as  that  of  the  wiggler  and  constant  values  of  | v^|  and 
, v„ |  are  allowed  by  the  system,  they  cannot  be  obtained 
with  arbitrary  inlet  conditions  in  the  electron  beam. 
Nevertheless,  these  are  the  trajectories  usually  used  in 


the  theory,4*8  without  studying  the  problem  of  how  the 
steady-state  situation  can  be  achieved.  An  additional 
complication  with  the  presence  of  an  axial  magnetic 
field  is  that,  as  will  be  shown  in  Sec.  II  of  this  paper, 
in  general,  there  exist  several  possible  steady- state 
trajectories  for  the  same  values  of  the  axial  field  and 
wiggler  parameters,  and  the  question  arises  as  to  which 
of  these  states  is  accessible  with  given  inlet  conditions 
of  the  electron  beam.  Thus,  in  the  presence  of  an  ax- 
iaL  guide  field  the  initial  conditions  and  the  structure 
of  the  transition  region  at  the  entrance  of  the  free  elec¬ 
tron  laser  may  be  of  crucial  importance  in  regard  to 
the  possible  modes  of  operation  of  the  device.  These 
factors  become  even  more  important  if  the  idea  of  re¬ 
cycling9  is  applied,  and  the  electrons  are  forced  to  pass 
the  transition  region  many  times. 

This  paper  presents  a  study  of  these  important  ques¬ 
tions.  In  Sec.  II,  we  derive  the  possible  steady  states 
in  the  homogeneous  part  of  a  free  electron  laser  and 
study  the  stability  of  these  states  to  perturbations  of 
electron  velocities.  The  transition  region  is  included 
in  the  theory  in  Sec.  Ill,  where  the  possible  ways  of 
operating  a  free  electron  laser  in  different  steady  states 
are  suggested  and  illustrated  by  numerical  examples. 

II,  EQUATIONS  OF  MOTION  AND  STABILITY 


Consider  a  cold  relativistic  electron  beam  moving  in 
a  magnetic  field  of  the  form 


B  =  B(z)el+  Vx  A  , 

a) 

where 

A  =-A(z)[excos<2>(z)  +  evsmd>(z)], 

(2) 

and 

d>  =  f  k0(z’)<W  . 

■'o 

(3) 

For  A  and  k0  independent  of  z,  the  vector  potential  (2) 
describes  the  field  on  the  axis  of  an  infinite  magnetic 
wiggler,  where,  as  is  well-known10 

A^I[pk0K0(pk0)  +  Kl(pk0)],  (4) 

where  /  is  the  current  in  the  wiggler,  p  is  its  radius, 
k0  =  2 ir/\t  A0  is  the  pitch  of  the  winding  of  the  wiggler, 
and  K0tl  are  the  modified  Bessel  functions  of  the  second 
kind.  By  using  the  more  general  form  (2)  for  the  vector 
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potential,  we  have  in  mind  primarily  the  possibility  of 
slow  variations  of  the  wiggler  parameters  p  and  k0  with 
z,  and  assume  that  in  this  case  the  magnitude  A(z)  in 
(2)  can  be  approximated  by  expression  (4),  where  p  and 
kQ  correspond  to  the  values  of  these  parameters  in  the 
nonuniform  wiggler  at  point  z. 


Although  the  magnetic  field  represented  by  the  poten¬ 
tial  (2)  does  not  satisfy  vxB  =  0,  it  gives  a  good  ap¬ 
proximation  of  the  exact  curl-free  field  on  an  infinite 
wiggler  at  small  distances  >*  from  it^  axis,  indeed,  as 
was  shown  in  Ref.  4,  the  relative  deviation  of  the  trans¬ 
verse  component  of  the  field  from  that  described  by  (2) 
is  of  the  order  of  (A’0r)2.  Accordingly,  if  the  beam  ra¬ 
dius  l  is  such  that  (k0l)2«  1,  the  actual  transverse  field 
can  be  well  represented  by  Eq.  (2).  The  axial  compon¬ 
ent  of  the  field  of  a  wiggler  near  the  axis  grows4  as 
feo r\  it  can  be  neglected,  however,  in  the  presence  of 
a  strong  axial  guide  field.  We  will  also  limit  ourselves 
to  low  current  beams  so  that  the  influence  of  the  self¬ 
space  charge  on  the  beam  can  be  neglected.  We  thus 
require  that  the  transverse  electrostatic  field  be  much 
smaller  than  vnBx/c  =  vnk0A/ct  or,  assuming  axial  sym¬ 
metry  of  the  beam,  w2pr«  2 eR0Av„/  me,  where  u)p  is  the 
plasma  frequency.  If,  for  example,  k  =  0.3  cm,  B  -  500 
G,  and  i>„  the  maximum  current  density  allowed  by 
the  model  will  be  approximately  250  A/cm2.  The  small 
signal  gain  in  a  free  electron  laser  at  these  conditions, 
however,  may  be  substantial,11  so  that  the  results  of  the 
present  work  could  be  important  in  current  experiments 


We  now  consider  the  momentum  equation  for  the  elec¬ 
trons  of  the  beam 


I(VV)=-^VXB- 


(5) 


On  using  normalized  velocities  ut  =vt/c  and  “time”  r 
=  cf  and  defining  £(z)  =eA/mc 2  and  Cl  -eB/  me2,  one  can 


write  (11)  in  the  form 

/.  n\ 

",  ="2(V'3 

(12) 

f.  flN  1  d£, 

u2=-Ul\k0u,--)+--u  „ 

(13) 

.  _*nl  1  dt, 

U, - - —  !(,  ■ 

3  y  1  y  dz  2 

(14) 

First  consider  the  homogeneous  case,  where  Cl  =  Cl0 
=  const,  &0  =  const,  and  £-£0=const.  In  this  case,  Eqs. 
(12)-(14)  have  a  particular  solution 

n1o=0,  uM-  const, 

and 


^  k0uw-Cl/y' 

which  with  the  conservation  of  energy 

i/V  =  i  -  i4o  -  *4 

defines  the  values  of  u70  and  u20.  The  question  arises  as 
to  how  this  steady-state  solution  can  be  achieved.  One 
can  answer  this  question  only  by  considering  the  transi¬ 
tion  region  of  the  wiggler,  where  £  and  k0  may  depend 
on  z.  Here,  one  would  expect  that  for  initial  conditions 
in  the  beam  ul  =u2  =0,  when  the  vector  potential  £  grows 
slowly  enough  with  z ,  the  velocities  u2  and  u3  would 
gradually  approach  their  steady-state  values  u2Q  and 
/<10,  ar.d  at  the  same  time  ul  remains  zero.  It  can  be 
shown,  however,  that,  in  general,  this  cannot  be  the 
case.  In  fact,  one  gets  from  (13)  that  if  iq(z)s0, 


(15) 

(16) 


where 

y  =  [1  -  (t’2/c2)]'1/2 . 


(6) 


.  1  rf? 

u*~yTzu'=u' 


d_ 

dz 


(17) 


Let 

ex(z)  =-er sin <t>  +eycos<i> , 

e2(z)  =-e,cos0  -  eysin0  ,  (7) 

®j(*i =e«  • 

Then,  the  cononical  model  vector  potential  is 

A  =  A(z)e2  ,  (8) 

and 

dA 

B=jB(2)e’_  °Aej-  (9) 

On  expressing  the  velocity  v  in  terms  of  the  orthogonal 
vectors  et,  e2,  and  e3  and  using 


and  therefore  «2  =  £/>,  which,  on  using  (12),  requires 
that  J2(z)s0.  Thus,  in  the  presence  of  an  axial  magnetic 
field  and  for  the  initial  conditions  on  u  considered  here, 
iq  cannot  remain  zero  in  the  transition  region.  The 
maximum  that  can  be  expected  is  that  the  component  iq 
in  the  transition  region  remains  small  in  comparison 
with  u2  and  m3.  When  this  is  the  case,  and,  in  addition, 
tq  remains  small  as  the  beam  propagates  in  the  homo¬ 
geneous  part  of  the  wiggler  we  define  the  beam  to  be 
stable  and  now  proceed  to  the  study  of  this  special  kind 
of  stability. 

First,  consider  the  homogeneous  region  of  the  device 
and  in  this  region  let 


de .  ,  rfe,, 

dt  =*°l'3C2’  dt 


one  can  rewrite  (5)  as 

=vVV'3- ~  (v2B +kavji) , 


dt  me 

dv,  e 

y^  =  -vV^+— 


y~dt~'7^(k°Vl^~V27^)' 


«i(r)  =  mt(t),  m2(t) -m20  +  w2(t)  ,  ?/3(t)  =hm  +  ?c3(t)  ,  (18) 

where  rq0  and  ?q0  are  given  by  (15)  and  (16),  and  ir/r) 
are  small  perturbations  to  the  steady-state  solution. 
Then,  on  linearizing  Eqs.  (12)—  (14)  one  gets 


uq  =  aiv2  +  bw3 , 

(19) 

(id 

w2  =  -  ate  l , 

(20) 

Wj  ^CtVy^  , 

(21) 

where 
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(22) 


<J  =  *oM30  “  “  *  ~3°  . 

0  30  y  y  «ao 

6=Vao-^fl=- 

°^  °  y  y  «30 

c=*0lo/y. 

Equation  (19)  then  gives 


(23) 

(24) 


wlz=atv2^bwi  =  -(a2  -  bc)wt .  (25) 

Thus,  the  necessary  condition  for  the  stability  of  the 
electron  beam  is 


a2  -  be  >0  , 


(26) 


or 


Further  study  of  the  stability  problem  must  involve  a 
knowledge  of  u20  and  Let  us  combine  Eqs.  (15)  and 
(16);  there  results 


1  _  ^2  / 1 Y  **3a  _ -  * 

30  \y/  (ux-to/yk0)*  r2' 

(28) 

This  equation  can  be  rewritten  in  the  form 

F-Fl  -  F2S0 , 

(29) 

where 

F 

;  u2  ’ 

(30) 

and 


1  + 


(Vr)2 


(Uy>-n/k0y)2 


(31) 


Assuming  U0/V)2<1-  1  /y2  for  y  large  enough,  the  func¬ 
tions  F'j  and  F2  have  the  general  form  shown  in  Fig.  1 
for  various  values  of  ft.  It  can  be  seen  from  Fig.  1(a) 
that  for  ft  =0,  there  are  two  solutions  for  corres¬ 
ponding  to  different  directions  of  propagation  of  the 


FIG.  1.  Schematic  of  the  functions  Fx  and  F2t  defining  various 
possible  steady-state  solutions  for  u  o- 


0.0  1.2  2.4  3.6 


4.8 


6  0  £l(rad/cm) 


FIG.  2.  The  real  positive  branches  of  hjq  vs  the  cyclotron 
freauenev  ft.  characterizing  the  iruide  magnetic  field. 


electron  beam.  In  the  presence  of  a  weak  axial  magne¬ 
tic  field,  there  exist  two  additional  solutions  for  u M  ad¬ 
jacent  to  the  resonance  velocity  w3r  =  ft/ kQy ,  as  shown  in 
Fig.  1(b).  If  one  continues  to  increase  ft,  a  situation  is 
reached,  where  again  there  remain  only  two  real  solu¬ 
tions  for  ui0  [Figs.  1(c),  (d)].  The  diagram,  where  all 
possible  real  positive  branches  of  ui0  are  presented  as 
a  function  of  ft,  is  given  in  Fig.  2  for  a  sample  case  in 
which  y  =  1.587,  k0  =  1.5708  cm-1,  and  £0  =0.3873. 

Let  us  now  find  the  frequency  ftcr  at  which  the  roots  2 
and  3  on  Fig.  2  become  complex.  This  transition  cor¬ 
responds  to  the  point  A  on  the  figure.  One  can  find  ftcr 
by  observing  that  the  function  F  in  Eq.  (29)  has  only  one 
bounded  maximum  at  the  point  «•,  such  that  Fv(:/‘)  =0, 
or 


ft 

*k0y’ 


where 


(32) 


It  is  now  clear  that  when  F  has  four  real  roots,  they 
are  contained  in  the  following  intervals:  [-1,0],  [0,«3J, 
[w3r, w* ],  [k*,+1],  On  the  ends  of  these  intervals  the 
function  F  changes  its  sign,  which  makes  it  easy  to  find 
the  four  roots  numerically.  It  is  also  clear  that  the 
roots  in  the  last  two  intervals  become  complex,  when 
i^n*)  =0.  Simple  algebra  then  leads  to  the  following  ex¬ 
pression  for  ftcr; 


Considering  our  sample  case  shown  in  Fig.  2,  Eq.  (34) 
gives  the  value  ftcr  =  G.763  rad  cm. 

We  now  return  to  the  question  of  stability.  It  is  clear 
that  inequality  (27)  (which  is  the  necessary  condition  for 
the  stability)  is  satisfied  for  branch  4  in  Fig.  2,  since, 
according  to  (15),  u2Q  is  negative  on  this  branch.  Simple 
analysis  also  shows  that  branch  2  is  stable,  since  the 
left-hand  side  of  (27)  on  this  branch  reaches  its  maxi¬ 
mum  value  of  1  only  at  ft  =ftcr;  branch  3,  in  contrast, 
is  always  unstable. 
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lit.  TRANSITION  REGION 


=u>x|0=0,  can  be  expressed  as 


The  Inequality  (27)  is  a  necessary  condition  for  sta¬ 
bility  of  the  electron  beam  in  the  homogeneous  part  of 
free  electron  laser.  This  condition  becomes  sufficient 
if  the  electron  beam  enters  the  homogeneous  region  with 
small  enough  component  ux  in  its  velocity.  We  now  pro¬ 
ceed  to  the  study  of  the  transition  region,  where  as  will 
be  shown,  special  experimental  steps  must  be  taken  in 
order  to  get  a  stable  electron  beam,  corresponding  to 
various  branches  on  the  diagram  on  Fig.  2.  Let  us  as¬ 
sume  that  the  vector  potential  £  in  the  transition  region 
is  a  slowly  growing  function  of  2.  Experimentally,  this 
would  be  the  case,  for  example,  if  one  gradually  de¬ 
creases  the  radius  p  of  the  wiggler,  or  increases  the 
pitch  length  =  2n,/&0  at  the  end  of  the  device  as  can  be 
seen  from  Eq.  (4).  Following  the  ideas  used  in  the  pre¬ 
vious  section,  one  can  find  approximate  solutions  of 
(12)-(14)  by  using  expansions  (18),  where  u20  and 
are  now  functions  of  z  and  correspond  to  the  components 
of  the  velocity  in  the  homogeneous  case  with  param¬ 
eters  such  as  those  at  the  point  z  in  the  transition  re¬ 
gion.  Then,  similar  to  (19)-(21)  one  has 


it,  =aw2  +  bw3 , 

(35) 

W2  *  -<HO,  +/«30  , 

(36) 

w}=cwl-fu20. 

(37) 

»x<* =77^)  siaMr>  -  HT')]dT'  •  (46) 

Thus,  if  the  vector  potential  £  grows  slowly  in  the  tran¬ 
sition  region  (the  function/  =  l/y  d£/dz  is  small  enough), 
one  expects  the  electron  beam  to  enter  the  homogeneous 
part  of  the  wiggler  with  a  small  magnitude  of  u\ ,  which 
is  sufficient  for  the  stability  of  the  beam  in  this  region 
if  the  inequality  (27)  is  satisfied. 

In  such  an  adiabatic  case,  one  can  also  find  the  tra¬ 
jectories  of  the  electrons  passing  the  transition  region. 
Expressing  the  radius  vector  r  of  the  electrons  in  the 
beam  in  terms  of  the  unit  vectors  el#  e2,  and  e3  [see 


Eq.  (7)],  one  has 

r  =  rlel  +  r2e2  +  ze3 ,  (47) 

which  on  differentiation  with  respect  to  r  gives 

r  =  (fl-  k^2r2)tl  +  (r2  +k0u2rl)e2  +  ze2 ,  (48) 

and  therefore  the  trajectories  are  described  by 

=Mi  +kou3r2  >  ^2  ~U2  *  *ow3ri  *  (49) 

This  system  of  equations  can  be  solved  in  the  following 
way:  Let 

R-rl+ir2i  Usul+iii2-  (50) 


where  a,  b ,  and  c  are  given  by  (22)- (24)  and 


(38) 


Then,  on  multiplying  the  second  equation  in  (49)  by  i  and 
adding  it  to  the  first  equation,  one  gets 

R  =  U-ik0u^R.  (51) 


Taking  the  time  derivative  of  Eq.  (35)  and  assuming  that 
the  coefficients  a  and  b  are  slGwly  varying  functions  of 
z ,  one  gets  the  following  equation 

=  aw2  +  btv2  +  aw2  +  bw3  =*  aiv2  +  bwz 

=  -(a*  -  bc)tvl  +f(au2Q  -  bu^)  ,  (39) 

or,  on  using  (22)- (24) 

=  -y?wx+g1  (40) 

where 

ii2  =  a2  -  be  ,  (41) 

and 


If  one  splits  R  into  two  parts 

R=R0+Rlt  (52) 

where 

iR„=U/k0u3,  (53) 

and  Rx  is  assumed  to  be  small,  then  on  linearization  in 

Eq.  (51), 

Rx  =-A0-  ik0u3Rl .  (54) 

The  solution  of  this  equation  for  Rx  is  given  by 

i?i=-  f  A0(T')exp(-i{(J)[£(r)]- <f>[2(T’)]})dT' ,  (55) 


*  =  (42) 

Assuming  that  p  is  a  slowly  varying  function  of  r  one 
can  approximate  the  solution  of  the  homogeneous  equa¬ 
tion 


it\  =-  P2k/  ,  (43) 

by  the  WKB  solution 

w,(r)  =  p'1'2[CiCos^(r)+Cssini('(r)],  (44) 

where 

■j)=  f  n(T')dr' .  (45) 

Jo 

Then,  it  can  be  easily  shown,  using  the  method  of  varia¬ 
tion  of  constants  in  (44),  that  the  solution  of  the  inhomo¬ 
geneous  equation  (40),  with  the  initial  conditions  w1|0 


where  0  is  defined  by  (3).  Thus,  if  the  velocities  ult 
u2t  and  Mj  are  slowly  varying  functions  in  the  transition 
region,  Rl  remains  small  along  the  trajectories  and 
R*R0,  or 


rx(z)^Re(R0) 


^ _ ffaoXsl — 

*o(*  )u^{z ) 


r2(z)^lm(R0)^ 


*0(Z)«30(Z) 


(56) 

(57) 


and  therefore  the  electrons  in  this  case  are  moving  on 
helical  orbits  with  adiabatically  changing  radius  rlf  and 
the  pitch  period  as  that  of  the  wiggler. 

Thus,  we  have  shown  that,  in  principle,  one  can  ob¬ 
tain  a  stable  electron  beam  in  a  free  electron  laser  if 
the  variations  of  the  parameters  of  the  wiggler  in  the 
transition  region  are  slow  enough.  This  conclusion, 
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&  =  0.74  rad/cm 


«  — 

/0.2  ■ 
/  0.2  ■ 
'  0.1  ■ 


20  30 


FIG.  3.  The  real  positive  branches  of  t*30  vs  4-  (a)  0 
=  0.6  rad/cm;  (b)  ft=  1.0  rad/cm,  (c)  ft  =  4.0  rad/cm. 


however,  is  based  on  the  approximate  solution  (46)  for 
wl  and  one  has  to  check  whether  all  the  assumptions, 
used  in  the  derivation  of  this  solution,  are  correct.  One 
of  these  assumptions  was  the  slowness  of  variation  of 
the  coefficients  a  and  b  in  (35)  as  the  beam  propagates 
through  the  transition  region.  Let  us  show  now  that,  in 
general,  this  is  not  guaranteed  even  if  4  varies  slowly. 
The  reason  is  that  the  real  solutions  for  and  w30, 
which  are  used  in  the  definitions  of  a  and  6,  do  not  al¬ 
ways  behave  continuously.  We  demonstrate  such  a  pos¬ 
sibility  in  Fig.  3.  In  this  figure,  one  can  see  the  dia¬ 
grams  of  the  possible  real  positive  solutions  for  uM9  ob¬ 
tained  in  a  fashion  similar  to  the  diagram  in  Fig.  2,  but 
for  constant  values  of  fi  and  varying  4.  Our  sample  case 
parameters  y  =  1.587  and  fe0  =  1.5708  were  again  used  in 
these  graphs.  As  mentioned  previously,  the  variation 
of  4  with  constant  value  of  kQ  can  be  experimentally  ob¬ 
tained  by  varying  the  radius  of  the  wiggler  winding  in 
the  transition  region,  holding  the  pitch  length  A0  =  2tt/£0, 
constant.  In  Figs.  3(a,b),  we  show  two  cases  with  the 
values  of  ft  higher  and  lower  than  the  critical  value  ftcr 
in  the  homogeneous  region  [ficr  is  defined  by  Eq.  (34), 
and  in  our  sample  case  is  equal  to  0.763  rad/cm].  For 
n<ftcr,  as  4  increases  in  the  transition  region,  one 
follows  the  path  AB  in  Fig.  3(a)  and  passes  continuously 
to  the  homogeneous  region  corresponding  to  the  point  B 
on  the  diagram  (at  this  point  4  =  40).  The  beam  is  stable 
in  this  case.  In  contrast,  if  ft  is  larger  than  ftcr,  one 
arrives  in  the  transition  region  at  the  point  D  [see  Fig. 
3(b)],  where  the  branches  2  and  3  of  u30  become  complex 
and  the  homogeneous  region  can  be  only  reached  on  the 
diagram  by  the  discontinuous  path  DEF.  The  jump  DE 
in  Ujq  leads  to  the  fast  variation  on  the  right-hand  side 
of  Eq.  (12),  which  cannot  remain  small  anymore,  and 
ux  grows  in  amplitude,  leading  to  the  instability  of  the 
beam.  For  sufficiently  large  values  of  0,  one  can  again 
return  to  the  stable  regime.  In  fact,  if  u2r*<l/yk0>  1, 


FIG.  4.  The  2  dependence  of  various  components  of  the  elecr 
tron  velocities  for  ft  =  0.74  rad/cm. 


the  only  possible  real  branch  of  is  branch  4  (see  Fig. 
2).  This  situation  is  illustrated  in  Fig.  3(c).  The  beam 
follows  a  continuous  path  GH  in  the  transition  region  in 
this  case  and  remains  stable. 

In  addition  to  these  qualitative  considerations,  we  il¬ 
lustrate  the  creation  of  the  instability  in  the  beam  in 
Figs.  4  and  5,  where  the  numerical  solutions  of  Eqs. 
(12)— (14)  for  ut  are  presented  for  our  sample  case  for 
two  values  of  ft  =0.74  and  0.77  rad/  cm  (recall  that  ft;r 
=  0.763  rad/cm).  We  assumed  in  these  calculations  the 
following  z  dependence  of  the  radius  p  of  the  wiggler 
winding  in  the  transition  region: 


(Po,  **0i 

P~  (  P0  +  (* .  *c°. 


where  p0=2.5  cm  and  zQ~  8  cm.  The  sudden  transition 
to  the  unstable  behavior  when  one  goes  from  Fig.  4  to 
Fig.  5,  where  all  the  ut's  change  rapidly  (w3  even  be¬ 
comes  negative  on  the  parts  of  the  trajectory)  is  ob¬ 
vious. 

Thus,  in  conclusion,  if  the  vector  potential  4  varies 


11  =  0.77  rad/cm 


FIG.  5.  The  z  dependence  of  wj  and  w3  for  ft  0.77  rad/cm. 
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slowly  in  the  transition  region,  one  can  get  a  stable 
electron  beam  for  ft<ftcr,  and  as  the  parameters  of  the 
wiggler  vary  adiabatically,  the  beam  follows  branch 
2  of  the  possible  solutions  for  on  Fig.  2.  One  can 
also  have  a  stable  situation  for  large  axial  magnetic 
fields,  when  branch  4  remains  the  only  possible  one  for 
operation.  One  has  to  remember,  however,  that  the 
necessary  condition  for  the  last  possibility  is  that  in 
the  transition  region  uZr  =  Sl/kQy  >1.  This  condition  can 
easily  be  satisfied  when  the  growth  of  £  in  the  transition 
region  is  due  to  the  variation  of  the  radius  p  of  the  wig- 
gier,  when  kQ- const.  If  in  contrast,  p=const  and  k0  is 
increasing  as  one  approaches  the  end  of  the  wiggler, 
larger  values  of  the  axial  magnetic  field  are  required  in 
order  to  operate  the  device  on  branch  4. 

Let  us  finally  consider  the  question  of  whether  it  is 
possible  with  the  initial  conditions  on  the  beam  assumed 
here  (namely,  ut  |0  =u2  j0  =0)  to  get  a  stable  electron 
beam  at  a  larger  region  of  branch  4,  especially  for 
a3r^  1.  As  mentioned  before,  the  necessary  condition 
(27)  for  stability  is  always  satisfied  on  this  branch, 
which  makes  it  more  attractive.  The  perpendicular 
component  of  the  velocity  on  branch  4  can  also  become 
very  large,  which  is  again  very  important  for  possible 
electromagnetic  wave  amplification  in  the  ^  direction. 

The  experimental  scheme,  which  allows  one  to  oper¬ 
ate  a  free  electron  laser  on  branch  4  is  shown  in  Fig.  6. 
We  are  exploiting  the  stability  of  the  beam  for  large 
values  of  ft  [as  demonstrated  in  Fig.  3(c)]  and  are  ap¬ 
plying  a  strong  axial  magnetic  field  in  the  transition  re¬ 
gion  of  the  wiggler.  Then,  after  passing  this  region  the 
electrons  will  enter  the  homogeneous  part  of  the  wig¬ 
gler,  being  on  the  upper  part  of  branch  4  in  Fig.  2.  Now 
in  the  homogeneous  region,  where  £  =  £0  =  const,  one  can 
gradually  reduce  the  axial  magnetic  field.  The  beam 
will  then  follow  the  continuous  branch  4  and  one  can 
easily  reach  region  ft=*ftcr,  which  was  unstable  with  the 
constant  axial  magnetic  field.  We  demonstrate  this 
possibility  in  Fig.  7,  where  the  solutions  of  Eqs.  (12)- 
(14)  are  shown  for  exactly  the  same  final  ft  and  £0  as  in 
the  unstable  case  in  Fig.  5.  The  same  variation  (58) 
for  p  was  used  in  the  computations.  The  cyclotron  fre- 


FIG.  6.  Possible  configuration  of  the  pump  and  guide  fields  for 
operating  on  branch  4  of  the  steady-state  regimes.  A— tran¬ 
sition  region  for  £;  B-transitton  region  for  ft;  C-homogenous 
part  of  the  device. 


FIG.  7.  The  z  dependence  of  the  electron  velocities  in  oper¬ 
ating  on  branch 4  with  varying  guide  magnetic  field.  0 
•  4  rad/cm,  0.77  rad/cm. 


quency  was  assumed  to  have  the  form 

L, 

°=  z>2L, 


(59) 


where  ft:  =  4  rad/cm,  ft2  =0.77  ra.d/cm}  and  L  =4  cm.  It 
can  be  seen  from  Fig.  7  that  the  beam  remains  stable 
and  corresponds  to  branch  4  with  negative  and  larger 
values  of  than  in  Fig.  4,  which  corresponds  to  branch 
2. 


IV.  CONCLUSIONS 

(i)  In  operating  a  free  electron  laser  with  an  axial 
magnetic  field,  different  steady-state  regimes  of  the 
helical  motion  of  the  electrons  in  the  homogeneous  part 
of  the  wiggler  must  be  considered. 

(ii)  The  necessary  condition  for  the  stability  of  these 
steady-state  regimes  is  given  by  the  inequality  (27). 

(iii)  The  transition  region  of  a  free  electron  laser 
plays  an  important  role  in  determining  the  sufficient 
conditions  for  stability  and  in  achieving  the  different 
modes  of  operation  of  a  free  electron  laser  for  a  given 
set  of  parameters  of  the  homogeneous  part  of  the  de¬ 
vice. 

(iv)  The  following  two  models  have  been  analyzed  for 
operating  a  free  electron  laser  in  different  steady- state 
regimes: 

(a)  The  first  is  characterized  by  a  constant  axial  mag¬ 
netic  field  and  gradual  increase  in  the  vector  potential 
in  the  transition  region.  The  stability  of  this  scheme  is 
limited  by  a  critical  value  of  the  axial  field  given  by 
Eq.  (34).  The  value  of  the  perpendicular  component  of 
the  velocity  is  also  limited  in  this  steady- state  regime. 

(b)  The  second  setup  uses  a  strong  axial  magnetic  field 
in  the  transition  region.  The  field  is  then  adiabatically 
decreased  in  the  homogeneous  part  of  the  wiggler.  This 
regime  is  always  stable  and  can  operate  with  any  value 
of  the  axial  magnetic  field  in  the  homogeneous  region. 
The  only  limitation  is  imposed  by  the  increasing  radius 
of  the  helical  trajectories  of  the  electrons  in  the  beam 
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as  the  perpendicular  component  of  the  velocities  grows 
with  a  decrease  in  the  axial  guide  field. 
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The  linearized  theory  of  a  free-electron-laser  amplifier  consisting  of  a  relativistic  electron  beam  transported  along 
the  axis  of  a  helical  wiggler  in  the  presence  of  an  axial  guide  field  is  solved  exactly.  With  suitable  re-idemification  of 
parameters,  the  theory  also  applies  to  the  case  where  the  wiggler  is  replaced  by  a  circularly  polarized  subiuminous 
radio-frequency  pump.  The  dispersion  relation  is  derived  and  numerical  examples  of  solutions  are  presented.  These 
indicate  (a!  that  the  use  of  an  axial  field  permits  operation  of  a  laser  of  given  high  frequency  and  undulatory 
transverse  velocity  of  the  unperturbed  electron  beam  at  lower  values  of  the  pump  field,  (b)  that  the  gain  can  be 
enhanced  by  approaching  the  condition  of  resonance  between  the  effective  frequency  of  the  pump  and  the  cyclotron 
frequency,  and  (c)  that  the  breadth  in  frequency  of  the  region  corresponding  to  spatially  exponentially  growing 
operation  can  be  much  extended. 


I.  INTRODUCTION 

The  theory  of  a  free-electron  laser  (FEL),  con¬ 
sisting  of  a  relativistic  electron  beam  transported 
along  the  axis  of  a  helical  pump  magnetic  field, 
has  been  given  by  Bernstein  and  Hirshfield.1  Their 
analysis  was  valid  for  arbitrary  pump  strength  but 
weak  rf  fields,  since  it  involved  linearization  in 
the  amplitudes  of  the  high-frequency  quantities. 
Here  we  present  the  extension  of  that  work  to  the 
case  where,  in  addition,  there  is  an  axial  magnetic 
field,  conventionally  present  for  beam  collimation. 
It  is  also  shown  that  with  a  suitable  reinterpreta¬ 
tion  of  parameters,  the  same  theory  applies  when 
the  magnetostatic  pump  is  replaced  by  a  circularly 
polarized  subiuminous  rf  pump.  The  axial  field  is 
shown  to  yield  the  additional  benefits  of  permitting 
the  use  of  weaker  pumps,  providing  enhanced  gain 
and  yielding  broader  domains  of  spacial  instability. 
This  is  discussed  in  detail  in  Sec.  VI. 

The  work  proceeds  as  follows.  The  general 
mathematical  description  is  developed  in  Sec,  n 
where  the  continuity  and  momentum  equations  de¬ 
scribing  the  relativistic  beam,  and  those  govern¬ 
ing  the  electromagnetic  fields  are  presented. 
Section  HI  describes  the  properties  of  a  helical 
pump  magnetostatic  field,  and  Sec.  IV  those  of  a 
circularly  polarized  subiuminous  rf  pump.  The 
linearized  equations  governing  the  high-frequency 
fields  are  derived  in  Sec.  V.  Section  VI  is  devoted 
to  a  brief  discussion  of  the  relation  of  this  work  to 
its  predecessors,  to  a  description  of  the  numerical 
examples  worked  out,  and  conclusions  concerning 
the  effects  of  the  axial  field. 


II.  GENERAL  MATHEMATICAL  DESCRIPTION 

Consider  a  cold  relativistic  electron  beam  de¬ 
scribed  by  the  continuity  equation 


gUv.(Nv)=0 

and  the  momentum  equation 


(nryv)  =  + 


**!)• 


(1) 


(2) 


where  m  is  the  rest  mass  of  the  electron,  and 

y=(l-r  Vc2)-1  2  (3) 

If  one  forms  the  scalar  product  of  (2)  with  7  v  and 
uses  (3)  to  express  v  in  terms  of  7  ,  there  results 
the  energy  equation 


r)/??0  =  -<*E  •  v 


(4) 


Let  £  be  a  constant.  It  is  convenient  to  introduce 
the  electromagnetic  potential  A  and  4>  via 

B  -Beg+  V  x  A  ,  15) 

z~vt- },(*)■  m 

Then  with  =eB/mc  one  can  write  (2)  in  the  form 
(±+v-?)(»/>v) 

=  +vx(V  x  A)  +  Bvx?,j 

smfiefxv+|^J  +  *  ^A-(rX)  •  ^  (7) 

or  on  rearranging  terms 

=  !Je.!<v  +  (— Wf4>_(VA)- v).  (8) 

*  \mc/ 

It  follows  from  the  Maxwell  equations 
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cVxB=4*J+^, 

(9) 

V.E*4t tZ  , 

(10) 

on  employing  (5)  and  (6),  that 

(11) 

^t  +  4itZ  =  -c'1  V.  — . 

at 

(12) 

Thus  if  we  adopt  the  canonical  model  of  FEL  the¬ 
ory,  viz. 


A  =A,U,f)iI+A,(2,f)e,  ,  (13) 

$=*(z,f)  (14) 


(note  that  the  vector  potential  is  written  in  the 
Coulomb  gauge)  and  assume  that  the  only  charged 
particles  present  are  electrons,  whence  2=-iVe 
and  J  =  -iVev,  then  (11)  and  (12)  yield 


(v- ete,- v)  , 


(15) 


a24> 

W 


=  4t rNe  . 


(16) 


III.  MAGNETOSTATIC  PUMP 

Consider  the  case  of  a  free-electron  laser  in 
which  the  pump  magnetostatic  field  is  generated 
by  helical  windings  and  the  self-fields  of  the  elec¬ 
tron  beam  are  negligible.  Then  in  cylindrical  co¬ 
ordinates  p,0,z  the  vacuum  magnetic  scalar  po¬ 
tential  x  will  be  helically  invariant,  viz. 

X=X(P,9-V).  (17) 

where  2ir/ft0  is  the  pitch,  and  will  satisfy  Laplace’s 
equation 


The  general  solution  of  (18),  regular  at  p==0,  on 
separation  of  variables  is  readily  shown  to  be 

x  =-Bz+  £  xJJ"lk0p)coa[m(t)  -  + 

(19) 

where  the  xm  and  are  constants  determined  by 
the  details  of  the  helical  windings.  Recall  that  the 
Bessel  function 


.»>-£ 


(s£)- 


s!  (m  +  s)! 


(20) 


Thus  if  a  is  the  radius  of  the  windings  and  p«  2ir/ 
fe0,  the  potential  is  well  approximated  by  the  term 
with  w*  1  alone,  with  Ix  approximated  by  the  lead¬ 
ing  term  in  the  series.  The  resulting  expression 


for  the  associated  magnetic  field  is,  on  choosing 
the  coordinate  system  so  that  =0  and  ^/?0x1 

B=-Vx-Be,+£0(e,cosfe0z  +  ey  sinAr^)  .  (21) 

The  nonconstant  part  of  (21)  can  be  written  as  the 
curl  of  the  vector  potential 

A  =  -(B0A0)(e,cosV  +  e,  sinA?^).  (22) 

Expression  (22),  valid  only  near  the  axis,  is  the 
form  conventionally  taken  for  the  magnetostatic 
pump  field.  A  corresponding  solution  for  the  ve¬ 
locity  and  density  can  be  obtained  from  (1)  and  (8) 
by  introducing  the  basis  vectors 


?!  =  -ex  sinAr 0z  +  ?,  cosA^  ,  (23 ) 

e2  =  -eIcos*az  -  e,sin/?oZ  ,  (24) 

e3=e,,  (25) 

when  on  writing 

A  =AieL+A2e2  +  A2e2  (26) 

it  follows  that 

3^  (*Al  t  ,  /QA2  ^  a\~ 

= -  kaAv 6i + (if + *«*•) e* 

+— ie3 . 

3z  3 

Thus  (1)  and  (8)  imply 

{h +v-£){yv~^)  -  k'”{yv*  -  e£) = > 

(it  +viP)  (yv>  ~  £) + k°v>{yv'  -  > 

(  e  \  (cat  a A.  ba.  .  .  .  .  \ 

=(  ~c  )  ("IT  ~Vy~at  ~V2Tf~  k°V2Al + k°v'At)  • 

(31) 

Now  on  combining  (22)  and  (25)  one  can  write 

A0  =  (mc2 /e)£0e2  ,  (32) 


(27) 

(28) 

(29) 

(30) 


where  £0  is  a  dimensionless  constant.  It  is  then 
readily  seen  that  if  also  4>o=0,  corresponding  to 
£o=-v*o=0,  then  a  solution  is  given  by 

v0  =  h  e3  +  we2  ,  (33) 

N0  =  const,  (34) 

where  u  =  const  and  =  const,  satisfy  (27)- (30) 
provided  that  consequent  to  (29) 

w  —  A?0  cu  £q(  A?gW  y  o  —  fl)  1  ;  (35) 
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where 

y0-[l  -(n^w*)/^]-4  2  (36) 

This  solution  and  its  experimental  accessibility 
has  been  analyzed  in  detail  by  Friedland.2 

IV.  RADIO-FREQUENCY  PUMP 

^The  solution  given  by  (32),  (33),  and  (34)  with 
Eo  =  0  can  also  be  adapted  to  describe  the  case  of 
a  free-electron  laser  with  an  electromagnetic 
pump  which  in  the  laboratory  frame  has  a  phase 
velocity  less  than  the  speed  of  light.  One  then 
views  the  solution  as  given  in  the  frame  where  the 
pump  wave  is  at  rest.  Equation  (15)  then  requires, 
on  using  (32)  and  (33),  that 

-k20c2Z0  =  J2pw/c ,  (37) 

which  on  using  (35)  can  be  written 

-k0c  =  J^ikguy 0  -  ft)'1  ,  (38) 


where 


Clearly  A'  is  a  Lorentz  invariant.  Also 
*V#  =V?(1  -  v0 u/c2) , 


(49) 


(50) 

(51) 

(52) 


w  = 


w/y 

1  -  VqU/C2  * 


(53) 


The  inverse  transformations  to  (50)- (53)  can  be 
gotten  by  interchanging  primed  and  unprimed  vari¬ 
ables  and  changing  the  sign  of  r0. 

The  counterpart  of  (35)  is  now 


.  (54) 

Equation  (38)  is  carried  into 


where  the  plasma  frequency,  defined  using  the 
rest  mass,  is 

u>,  =  (4t iNQe?/mY2 .  (39) 

Let  vQ  be  the  speed  of  the  laboratory  frame  as 
seen  from  the  wave  frame.  Distinguish  quantities 
in  the  laboratory  frame  by  a  prime.  Then  on  Lo¬ 
rentz  transformation  2'  =  >(c  -  r0t)9  t*  =  >  (/  -  v0z:' 
c2),  and 

w'  =-k0v0y  ,  (40) 

k'0  =  Ky  >  (41) 

where 

f-U-^/c*)-1'*1  .  (42) 

Clearly 

v0  =_w'/*4,  (43) 

is  the  negative  of  the  phase  velocity  of  the  wave. 
Moreover, 

BosE^  +  ^otexCOS^'-  w't')  +  ey  sin(A£z' -  u>'f')] , 

(44) 

where 

&o  =  yB0  (45) 

and 

E'  =  -(oj'/^c)e3xg'.  (46) 

Evidently  the  wave  is  transverse  and  circularly 
polarized.  The  associated  potentials  are 

*o  =  0,  (47) 

Ao  =A'[et cos(fc£z'  -  uj'/')  + e,  sln(k^z’  -  w'/')] , 

(48) 


-  fe'V)  =  ^(^/'  -  u>)[y’0(k'nu>  -  u.-')  -  ft]'1  . 

(55) 

Equation  (55)  can  be  viewed  as  the  dispersion  re¬ 
lation  for  the  pump  electromagnetic  field,  but  it 
is  to  be  noted  that  the  steady- state  theory  is  not 
restricted  to  weak  pump  fields  and  a  linearized 
theory. 

V.  STABILITY  ANALYSIS 

Let  us  work  in  the  laboratory  frame  for  the  case 
of  the  magnetostatic  pump  and  in  the  wave  frame 
for  the  case  of  the  radio- frequency  pump.  The 
stability  analysis  is  then  common.  Let 

A  =  A„  -  Re{( m c* /e)[ S, U )e, U )  +  $2(z  )e2(£ )  ]<?'*" »}  , 


(56) 

vsVo  +  RetVUJc'1"*] ,  (57) 

4>  =  0  +  R  e[(mcVe)(w/*c)5,e-,"‘] ,  (58) 

y  =  y0  +  Re(r  e~iu‘) ,  (59) 

N  =N0  +  Re(.Vle~iwl) .  (60) 

Then  (29)  and  (30)  yield  on  linearization 

(_“*H‘“5«)(r  0V,  +  C^)-*0“(rM,  +  >0^2+C^) 

-k0V,{y0u-c£,0)  =  -QV2  (61) 


(-.■w+g)(rH+>0v2+cy+v(>°v71  +  ^i)=nv'l . 

Rather  than  use  (31)  it  is  convenient  to  employ  the 
linearized  version  of  (4)  which  yields 


(62) 
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Linearization  of  (1)  gives 


(-'“*sMs)(w=°. 

(63) 

while  linearization  of  (3)  implies 

r/yl*(uV3+wV2)/c*. 

(64) 

Equations  (15)  and  (16)  on  using  (27)  yield 

^•-^-2*0^+^=--^, 

(65) 

(66) 

kc  dzi  Ir  N0  * 

(67) 

Note  that  Eqs.  (60)-(67)  are  a  system  of  eight  lin¬ 
ear  ordinary  differential  equations  with  constant 
coefficients  for  the  eight  quantities  Vlt  V2f  K3)  r, 
Nlf  £lt  £2,  and  £3.  Thus  we  may  seek  solutions 
where  all  these  scalars  vary  with  z  as  eik *.  If  we 
write 

E=Re[aU)ei“‘],  (68) 

then  it  follows  from  (6)  that 

a  =  -i(mcw/e)| .  (69) 

Equations  (60)  through  (67)  then  imply 

i(kll  -  ^)(>  o^l  +  c?l)  -  -  ^/>'o)(Vo^2+  C^2+  TlV) 

=  (fi/y0)c?2+  (to/y0)rw+  k0V3(y0w  -  cy  ,  (70) 

(kj4-tt/y0)(y0V1  +  cti) 

+  i(ku-  w)(y0V2  +  c?2+  Tu;)  =  -(«/y0)c|l ,  (71) 

r=-  u*i+£*t ,  (72) 

c  ku  -  w 


kV3 

lV0  "  -  U>  ' 

r  =  (yJ/c2)(MV3  +  MJ^)f 

[1  -  c»(^+  **)/«•&  -  (2ik0k<?/ <i^)?2 

--(mJVXVc)  . 

(2»fe0*ci/«2)?l  +  [1  -  c2(J^+  *W]?2 
=  -(t4/uJ2)[(V2/c)  +  (it)/c)(N1/N0) , 
.V,AV0  =  _(*cw/c4)«3. 


(73) 

(74) 

(75) 

(76) 

(77) 


It  is  convenient  to  express  r,  iV lf  V\,  V2,  and  V*3 
in  terms  of  {lf  $2,  and  ?3.  The  result  can  be  rep¬ 
resented  in  the  form 


£•1=0, 

where  the  dielectric  tensor 


(78) 


<_=£+  (u)*pT  /VoW2)* 


(79) 


r  =(n/r0)[(if«-  u>)2-  (V-  n/>0)2]-1 .  (80) 

The  components  of  0  are 

0U  *  1  -  c2(k^,+  &)/J*  -  wj/v0^2  , 

e2J  =  i-c2(*2+fe2)/u^  (8i) 

-(^/Vo^K1  +  -  (I))-2] , 

023=8Jt  =  -(Jtt/yo)(ku-  w)**(w/w)(fe-  wk/c2)  , 
eM  =  i  -  (u>*/Vo)(*«  -  wj-'d  -  uVc2) . 


The  elements  of  tfi  are 

|Jlu  =  feou-%'o . 

kt;). 

— s — "/> 

.  /  fc^  +  w2)  ,  u ?  U?  \ 

- u - w/  ku-*  X L  +  ?r  ’ 

/fc(K2  +  W2)  \  /?0M  -  Q  hn 

^23 "  \  u  /  ku-  U) 


(82) 


/ttW  <«>  \ 

\"?r  ku  -  u)  c4  -  ft/y0  /  ’ 


#31  =  iww/u  , 


uw  -  fl/y0 
M  ku - 


ku-u/' 


_u rw  k0u-tt/yQ/uw  u)  Z0k0cw  ku  -  u  \ 

^33  “  m  ku  -  w  y'c3’  -  u>  +  wj  /?qU  -  Q/y0/ 

In  the  limit  12-0,  r  vanishes  and  £  reduces  to  £ , 
which  apart  from  notation  is  the  form  found  by 
Bernstein  and  Hirshfield.1 


VL  THE  DISPERSION  RELATION  AND  NUMERICAL 
EXAMPLES 

In  order  that  (78)  have  nontrivial  solutions  it  is 
necessary  that  the  determinant 

D  “ det£  =  0  .  (83) 

This  yields  an  eighth-order  polynomial  equation 
for  k.  In  practice,  for  the  cases  of  interest  u£ 

«  u j2  and  m~c,  and  two  of  the  roots  aie  such  that 
uj/fes-c.  That  is,  they  propagate  in  the  negative- 
z  direction  counter  to  the  beam  and  are  substanti¬ 
ally  unaffected  by  the  tenuous  beam.  The  remain- 
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ing  six  roots  correspond  to  waves  which  propagate 
along  the  beam.  When  the  two  of  these  which 
can  be  associated  with  cyclotron  waves  in  the  limit 
of  no  helical  pump  disappear,  and  one  recovers 
the  result  of  Bernstein  and  Hirshfield.1  These  fea¬ 
tures  will  be  illustrated  later  when  numerical  ex¬ 
amples  are  discussed. 

Now  Eqs.  (61 )— (67)  comprise  a  tenth-order  sys¬ 
tem  of  linear  ordinary  differential  equations  which 
require  for  a  unique  solution  the  stipulation  of  ten 
boundary  conditions.  Since  usually  there  is  negli¬ 
gible  reflection  of  waves  at  the  output  end  of  an 
FEL  amplifier  of  finite  length,  two  conditions  are 
the  requirement  that  the  amplitudes  of  the  waves 
propagating  counter  to  the  beam  be  zero.  This 
requirement  can  be  most  easily  dealt  with  via 
solving  the  system  of  ordinary  differential  equa¬ 
tions  by  means  of  a  Laplace  transform  in  z,  as 
was  done  in  Ref.  1,  instead  of  the  normal  mode 
analysis.  The  dispersion  relation,  of  course,  de¬ 
termines  the  poles  of  the  transform  in  terms  of 
which  the  inversion  can  be  readily  accomplished. 
The  resulting  solution  for  a(z)  can  be  written  in 
terms  a(0),  assuming  that  all  other  first-order 
quantities  are  zero  at  z  =0  and  involve  linear  com¬ 
binations  of  the  six  modes  corresponding  to  the 
six  roots  with  Refc>0.  Since  in  general  these  roots 
are  nondegenerate,  but  differ  by  amounts  of  order 
much  less  than  x’V,  there  will  be  interference 
amongst  their  contributions  to  aU),  which  be¬ 
comes  evident  after  a  distance  of  order  2 
This  feature  has  been  examined  in  detail  in  Ref. 

1.  We  will  not  pursue  it  further  here,  other  than 
to  note  that  the  single  particle  theory  in  which  one 
examines  the  second -order  energy  change  in  a 
distance  z  of  an  electron  moving  in  the  zero-  and 
first-order  electromagnetic  field,  and  identifies 
this  with  the  gain  in  energy  of  the  high  frequency 
field,  is  valid  only  for  z$.k<  1. 

We  now  consider  the  dispersion  relation  (83)  in 
an  FEL  with  guide  magnetic  field.  Because  of  the 
complexity  of  the  dielectric  tensor  €  [see  Eq.  (79)] 
it  is  convenient  to  study  the  dispersion  relation  by 
comparing  two  FEL’s,  identical  except  that  one 
has  an  axial  field  while  the  second  does  not  and 
thus  is  characterized  by  the  dispersion  relation 
X)0  =det(8)-0,  the  properties  of  which  are  well 
understood.  We  make  the  comparison  between  the 
two  lasers  by  fixing  the  parameters  of  the  FEL 
without  the  guide  field  and  adjusting  the  value  of 
the  pump  field  parameter  $0  in  the  laser  with  the 
guide  field  so  that  the  axial  velocities  u  (and  there¬ 
fore  also  w)  in  both  lasers  are  identical.  This 
assures  the  same  Doppler  upshift  of  the  frequen¬ 
cies  in  the  lasers.  A  similar  comparison  has  been 
made  by  Friedland  and  Hirshfield  for  the  single 
particle  model  of  FEL.3 


Let  be  the  pump  field  parameter  in  the  FEL 
without  the  guide  field.  The  unperturbed  electron 
velocity  components  are  then  given  by  w/c  =  ?o/':u 
and  u/c  =  [l  -  (1  +  Q)/rl]1  2.  Therefore,  following 
Eq.  (35),  with  the  guide  field 


This  equation  demonstrates  the  intriguing  possibil¬ 
ity  of  reduction  of  the  pump  field  in  an  FEL  as  one 
approaches  the  cyclotron  resonance  condition  8/>0 
—  kf/i .  Accessibility  of  the  resonance,  however, 
is  not  guaranteed,  as  was  shown  in  the  recent 
study2  of  the  unperturbed  electron  beam  orbits  in 
an  FEL  with  the  guide  field.  It  was  demonstrated 
that  for  given  values  of  y0,  k0 ,  40,  and  ft  the  elec¬ 
trons  can  possess  more  than  one  steady  state. 

For  example,  Fig.  1  shows  u/c  versus  D/c  for 
fc0-6  cm*1,  y0  =  3,  and  £o  =  0.5.  For  it  is 

seen  that  only  one  branch  exists  (branch  C).  But 
when  D<ftcr  two  additional  branches  (A  and  B)  are 
allowed.  It  was  also  shown  that  the  necessary  con¬ 
dition  for  orbital  stability  of  the  steady- state  solu¬ 
tions  against  small  perturbations  is  given  by  the 
inequality 


Branch  C  is  always  stable,  since  /r<0  on  this 
branch.  On  branches  A  and  B ,  /r>0,  but.  as  was 
shown,  only  branch^  satisfied  (85)  and  thus  may 
be  used  in  applications.  Since  the  ratio  tv/u  is 
kept  constant  in  our  comparative  study,  one  can 
substitute  the  expression  for  $0  found  from  (35) 
into  (85)  and  write  the  stability  condition  in  the 
following  form: 


tl<Qrr  = 


1  +  ( iv  /u)2  1 


valid  for  branches  A  and  B .  In  our  sample  case 
(y0-3,  “ 6  cm*1,  and  ^  =  0.5)  one  has  S7cr/c 

=  16.18  cm*1,  and,  therefore,  according  to  (84), 


0  4  8  12  16  20  24  28 


fl/e  (cm*') 

FIG.  1.  Steady-state  normalized  axial  velocity  u/c 
as  a  function  of  normalized  axial  magnetic  field  n/c. 
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FIG.  2.  Dispersion  function  D*  on  branch  A  for  the 
case  y0  *3,  *0  =  6  cm*1,  and  u)/c  =40  cm*1.  The  dashed 
curve  represents  the  FEL  without  the  guide  magnetic 
field  and  £§=*  0.5.  The  solid  curve  is  for  the  FEL  with 
the  guide  field  (0/c=6.5  cm-1),  where  smaller  values 
of  £§  are  used  so  as  to  provide  the  same  values  of  u  and 
w  as  for  the  dashed  curve. 

£0  on  branch  A  cannot  become  less  than  ?cr  =  1.562 

xicr2. 

We  return  now  to  the  study  of  the  dispersion  re¬ 
lation  (83).  The  form  of  the  dielectric  tensor  c 
[Eq.  (79)]  suggests  that  for  values  of  u)p  small 
enough,  the  function  D  will  differ  significantly 
from  D0  only  in  the  regions  where  (ku  -u>)2  -  (fe0 u 
-tt/y0)2  31  as  a  result  of  the  resonance  in  the 
denominator  in  r  [see  Eq.  (80)].  We  demonstrate 
a  typical  effect  of  the  axial  guide  field  on  the  dis¬ 
persion  function!)  in  Fig.  2,  where  the  function 
D'=D[ku  -  u))/up{l  -  u2/c2)]2/y0  (the  full  line)  is 
shown  versus  n  =  ck/u>  for  branch  A  in  the  sample 
case  when  w/c  =  40  cm’1,  u^/c^O.S  cm"2,  and 
0/c  =  6.5  cm’1.  In  the  same  figure  the  dashed  line 
represents  the  case  with  no  guide  field. 

It  is  well  known1  that  the  unstable  regime  in  an 
FEL  without  the  guide  field  can  be  described  as  a 
coupling  between  the  transverse  electromagnetic 
modes  with  the  dispersion  relation  nx  a  =  1  ±  ckQ/(jj 
and  the  electrostatic  beam  modes  characterized 
by  n3  4  =  c/m  ±  cuj#/V0u«<.  One  can  see  from  Fig.  2 
that  these  four  roots  are  only  slightly  perturbed 
by  the  presence  of  the  axial  field.  There  exist, 
however,  two  additional  roots  in  the  neighborhood 
of  the  resonance  points  6  -c/u  ±  (c/?0/w  -ftc/ 
yuv).  If  the  resonances  are  widely  separated  as 
in  the  case  of  Fig.  2,  the  onset  of  the  unstable 
mode  is  roughly  the  same  as  without  the  guide 
field,  namely,  as  the  frequency  u>  increases,  the 
root  nx  moves  to  1,  passing  the  region  nA<n<rui 
(since  n3  4-c/w).  The  modes  couple  in  this  region, 


FIG.  3.  Graphical  representation  of  the  dispersion 
function  on  branch  A  for  the  case  y0  =  3,  k  0  =  6  cm-1, 

|J  =  0.5,  y/cs50cm“1,  and  Increasing  values  of  the 
guide  field  (the  solid  curves):  (a)  S7/c  =  14  cm”1,  (b) 
Q/c  ~  15  cm"1 ,  (c)  ft/c  =  16  cm-1.  The  dashed  curves 
correspond  to  the  FEL  without  the  guide  field.  Two 
pairs  of  roots  of  the  dispersion  relation  become  com¬ 
plex  as  the  real  roots  n3  and  n*  are  squeezed  by  the 
resonances  at  n$  and  n$. 


and  the  roots  of  the  dispersion  relation  are  com¬ 
plex.  When  continues  to  increase,  nx  becomes 
less  than  n4,  the  coupling  diminishes,  and  one 
again  has  a  stable  regime. 

New  effects  may  occur  when  the  resonances  6 
approach  each  other.  This  situation  is  shown  in 
Fig.  3,  where  the  full  line  represents  the  disper¬ 
sion  function  on  branch  A  for  increasing  values  of 
8.  One  can  see  in  this  example  that  even  for  <jl?/c 
=  50  cm*1  in  our  sample  case  (all  the  modes  are 
stable  in  this  case  if  8=0)  it  is  possible  just  by 
changing  8  to  squeeze  the  roots  w3  4  by  the  reso¬ 
nances  n5  6  so  that  two  pairs  of  the  roots  become 
complex.  For  higher  frequencies,  when  again  the 
FEL  without  the  guide  field  is  stable  (KiOi4)  one 
can  also  get  an  unstable  regime  as  is  demonstrated 
in  Fig.  4  for  uj/c  =  100  cm*1.  Our  numerical  study 
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FIG.  5.  Spacial  growth  rates  Im(fc)  versus  u/c  on 
branch  A  (solid  curves)  and  C  (dashed  curves)  for  vari¬ 
ous  values  of  Q/c :  (l)ft/c=0,  (2)  U/c  =12  cm"1 ,  (3) 
ft/c  =  14  cm*1,  (4)  fi/c  =  14.5  cm*1 ,  (5)  ft/c  =28  cm'1 , 
(6)  ft/c  =  23  cm"1,  (7)  l/c~  21  cm"1.  For  all  the  cases 
7o  *3,  =6  cm"1,  and  ^  =0.5. 


FIG.  4.  Graphical  representation  of  the  dispersion 
function  on  branch  A  for  the  case  y0  =  3,  k  o"6  cm"1, 
cj  =0.5,  o>/c  =100  cm"1.  The  solid  curves:  (a)fi/c  =  14 
cm’1;  (b)  ft/c  =  15  cm"1 ;  (c)  ft/c  =16  cm"1 ,  The  dashed 
curves  correspond  to  the  FEL  without  the  guide  field. 

shows  that  similar  behavior  is  also  characteristic 
for  branch  C  with  the  only  difference  that  there  is 
only  one  pair  of  unstable  modes  in  the  low  and  the 
high  frequency  ranges,  respectively. 

We  finally  summarize  our  comparison  of  the 
FEL’s  with  and  without  the  guide  field  in  Figs.  5 
and  6,  where  the  imaginary  part  of  k  is  shown  as 
a  function  of  u )/c  for  various  values  of  the  axial 
field  in  our  sample  case  (y0  =  3,  =  6  cm’1,  ^ 

=  0.5,  ujJ/c^O.S  cm2).  Figure  5  is  for  0 <ft/c 
<14.5  cm*1  on  branch  A  (the  full  lines)  and  21<ft/c 
<28  cm*1  on  branch  C  (the  dashed  lines).  The 
resonances  n5t6  are  relatively  wide  apart  from 
each  other  and  formally  the  instability  in  this 
range  of  ft  occurs  similarly  to  the  case  of  the  las¬ 
er  without  the  guide  field.  Nevertheless,  the 
presence  of  the  guide  field  increases  the  instability 
on  branch  A  and  tends  to  decrease  it  on  branch  C. 
In  addition,  the  linewidth  of  the  unstable  regime  is 
seen  to  be  significantly  increased  at  lower  fre¬ 
quencies  on  branch  A.  Together  with  this,  no  in¬ 
stability  exists  at  frequencies  higher  than  those 


characteristic  of  the  FEL  without  the  guide  field. 
As  one  approaches  the  resonance  condition  ft 
=y0uk0  (further  increasing  ft  on  branch  A  or  de¬ 
creasing  it  on  branch  C)  a  completely  different 
type  of  behavior  is  observed  as  is  shown  in  Fig.  6 
for  ft/c  =  15.25  cm'1  on  branch  A  (the  full  line)  and 
ft/c  =  18  cm’1  on  branch  C  (the  dashed  line).  The 
unstable  region  extends  over  the  entire  low-fre¬ 
quency  range  and  there  are  two  different  unstable 
modes  on  branch  A,  as  was  mentioned  previously. 
In  addition  there  exist  unstable  modes  in  the  high- 
frequency  region,  which  was  totally  stable  before. 
Note  that  the  values  of  Imfc  in  this  high-frequency 
regime  are  only  weakly  dependent  on  the  frequency 
itself. 

Thus,  in  conclusion,  we  have  demonstrated  in 


FIG.  6.  Spatial  growth  rates  Im(Jt)  in  the  sample  case 
(y0=3,  fe0=6  cm-1,  ^=0.5)  versus  ^/c  in  the  regime, 
where  the  cyclotron  modes  couple  to  the  beam  modes 
(see  Figs.  3,4).  Branch  A  (solid  curves):  Q/c  =  15.25 
cm*1.  Branch  C  (dashed  curves):  fi/c*18  cm-1.  The 
unstable  modes  are  extended  over  the  low-  and  high- 
frequency  regions.  There  exist  two  different  growth 
constants  in  this  regime  on  branch  A. 
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our  numerical  examples  that  the  presence  of  the 
guide  field  in  an  FEL  introduces  the  following  de¬ 
sirable  features: 

(i)  One  can  operate  the  laser  with  much  lower 
magnitudes  of  the  pump  field  without  sacrificing 
the  undulatory  velocity  of  the  electrons.  This  al¬ 
lows  one  to  use  shorter  periods  of  the  wiggler  with 
the  same  currents. 

(ii)  The  laser  can  be  operated  in  higher  gain  re¬ 
gime  by  approaching  the  resonance  condition  ft 


(iii)  The  linewidth  of  the  unstable  modes  can  be 
widely  extended  to  both  low-  and  high-frequency 
ranges. 
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Helical  magnetic  wigglers  for  free  electron  lasers  can  pro¬ 
duce  non-helical  electron  trajectories  if  a  uniform  axial 
guide  magnetic  field  is  imposed.  F reidland ’ 3  necessary 
criterion  for  the  existance  of  helical  orbit's  is  reviewed 
and  shown  to  apply  for  non-relativis tic  electron  energies. 
An  experiment  designed  to  test  this  criterion  is  described 
and  results  are  compared  with  theory. 

Key  words:  free  electron  laser,  magnetic  wiggler,  elec¬ 
tron  orbits. 


Introduction 


Considerable  effort  is  currently  underway  in  the  anal¬ 
ysis  (1),  design  (2),  and  construction  (3)  of  free  elec¬ 
tron  lasers  for  amplification  of  infrared  and  far  infra¬ 
red  radiation,  A  typical  device  comprises  a  good  quality 
electron  beam  with  energy  of  10 ?s  of  MeV  which  moves 
through  a  periodic  static  pump  magnetic  field-,  termed  a 
magnetic  wiggler.  Radiation  propagating  along  the  elec¬ 
tron  beam  has  been  shown  experimentally  (4)  to  be  ampli¬ 
fied,  but  the  single-pass  small-signal  gain  may  be  quite 
small  (7%  was  reported  for  a  320  cm  length  at  X  =  10. 6u 
in  Ref .  4) . 
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Suggestions  for  enhancing  the  small-signal  gain  by 
superposing  a  uniform  axial  magnetic  upon  the  wiggler  field 
have  appeared,  based  upon  both  single-particle  (5,6)  and 
collective  (7)  models.  The  gain  enhancement  can  result 
from  either  increased  equilibrium  undulatory  momentum  (5) , 
or  from  dynamical  resonance  between  induced  electromagnetic 
perturbations  and  the  natural  oscillations  of  electrons  on 
helical  orbits  (6,7).  The  increased  undulatory  momentum 
results  in  a  decreased  axial  momentum,  and  thus  a  decreased 
Doppler  up-shift,  i.e.  the  laser  output  frequency  is  shif¬ 
ted  to  longer  wavelength.  Gain  enhancement  may  still  be 
achieved  without  this  wavelength  increase  by  operating  the 
device  with  a  reduced  wiggler  field. 

A  necessary  condition  for  achievement  of  the  gain  en¬ 
hancement  is  that  the  equilibrium  electron  orbits  -in  the 
wiggler  be  nearly  helical.  Without  the  axial  guide  field 
a  helical  magnetic  wiggler  produces  a  helical  orbit;  this 
result  follows  from  the  constancy  of  canonical  angular 
momentum.  But  when  the  axial  guide  field  is  present,  the 
orbits  are  generally  not  helical  (8).  They  can  be  arranged 
to  be  nearly  helical  if  the  entry  conditions  into  the 
wiggler  are  suitably  tailored,  and  if  the  wiggler  and  guide 
field  parameters  are  in  a  regime  of  stability,  determined 
from  the  orbit  parameters  (9). 

In  this  paper,  we  shall  review  the  basis  underlying 
the  criterion  for  orbit  stability,  and  shall  present 
results  of  an  experiment  designed  to  test  this  criterion 
quantitatively. 


Orbit  Stability 


Here  we  summarize  (8)  some  aspects  of  the  dynamics  of 
charged  particles  moving  in  a  static  magnetic  field  given 
by 


B(z)  B  +  (£  cosk  z  +  &  sink  z)B, 
—  '  zo  x  o  •  y  o  -L 

=  §~B  -  e«B 

3  o  2 


(1) 


Here  B0  is  the  magnitude  of  the  uniform  axial  guide  field, 
and  B-l.  is  the  magnitude  of  the  transverse  helical  field 
with  pitch  £Q  *  2ir/k0.  It  has  been  shown  (1)  that  the 
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charged  particle  dynamics  in  this  field  are  described  com¬ 
pactly  if  a  coordinate  system  with  basis  vectors  (£]_>  ^2> 
£3)  is  used,  rather  than  the  Cartesian  system  (£x,  ev,  e2). 
The  coordinate  transformations  follow  from  the  definitions 
of  §2  anc*  ^3  given  in  Eq.  (1),  and  by  =  e 2  *  £3. 

Of  course,  the  field  given  by  Eq.  (1)  does  not  satisfy 
V  x  b  »  0;  it  is  however  a  good  approximation  to  the  actual 
field. near  the  axis  of  two  identical  interspersed  helical 
conductors  carrying  currents  in  opposite  directions.  The 
exact  field,  and  the  precise  nature  of  the  approximations 
leading  to  Eq.  (1)  will  be  discussed  in'  a  forthcoming 
paper  (10) . 

For  a  particle  of  charge  e,  rest  mass  m,  and  relativ¬ 
istic  energy  factor  y,  the  steady-state  solutions  of  the 
equation  of  motion  md(yyO/dt=  -ev^  x  1$  [with  _B  given  by 
Eq.  (1)]  are 


u  =  0 


= 


= 


ko5u3 


k  u^y  -  ft/c 


o  3 

2  -2,1/2 

(1  -  u2  -  Y  ) 


(2) 


where  u.  =  v/c,  ft  =  eB0/m,  and  £  -  eBj-/k0mc.  These  compo¬ 
nents  correspond  to  ideal  helical  trajectories,  since  U2 
and  U3  are  constants.  However,  these  steady-state  values 
can  only  be  approached  asymptotically,  for  an  actual  wig- 
gler,  because  of  coupling  between  the  components  in  the 
transition  region  at  the  entrance  to  the  wiggler  (8) ,  and 
because  the  form  given  by  Eq.  (1)  is  only  an  approximation. 


The  solutions  given  by  Eq.  (2)  are  depicted  (for  y  * 
10.0,  kQ  =  6.0  cnfl,  and  £  =  1.0)  in  Fig.  1.  For  ft  >  Qcr 
the  equations  are  single-valued,  whilst  for  ft  <  ftcr  they 
are  triple-valued.  The  critical  axial  guide  field  cyclo¬ 
tron  frequency  ftcr  is  given  by 


ft 

cr 


=  k.Qc[(Y2 


,,1/3  ,2/3, 3/2 

l;  ~  c,  j 


(3) 
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Figure  1.  Solutions  for  steady-state  axial 
momentum  U3,  as  a  function  of  axial  magnetic 
field,  for  y  »  10,  k0  =  6  cm“^,  and  5=1* 


For  smaller  values  of  5  than  that  chosen  for  Fig,  1,  the 
curves  hug  more  closely  the  asymptotes  U3  =  (1  -  y-2)l/2 
and  U3  =  ft/k0cy.  Perturbation  theory  shows  (8)  that 
branches  A  and  C  in  Fig,  1  are  stable,  whilst  branch  B  is 
unstable.  Thus,  if  a  particle  enters  a  wiggler  along  a 
gradually  increasing  guide  field,  it  would  move  on  a 
stable  helical  orbit  along  branch  A,  but  at  ft  =  ftcr  the 
orbit  would  become  unstable  and  thus  severely  non-helical. 
Examples  of  non-helical  orbits  are  shown  in  Ref.  8.  If 
ft  *  const  and  the  wiggler  field  increases  gradually,  a 
similar  phenomenon  occurs  at  £cr,  where 


(«  -  Kv2  -  »1/3  -  w/k  c)2/V/2 


or  (B/B)  =  [ Cy2  -  l)1/3(k  c/fl)2/3  -  1]3/2 

^  o  cr  o 


(4) 

(5) 
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Thus  for  a  charged  particle  moving  through  a  wiggler 
in  a  uniform  axial  guide  field,  the  orbit  can  be  nearly 
helical  if  5  <  £cr  a H  along  the  wiggler  but  would  depart 
significantly  from  helicity  if  g  >  Ccr* 


Experiment 


Although  electron  beams  of  interest  for  practical  free 
electron  lasers  have  relativistic  energies,  the  phenomenon 
of  helical  orbit  stability  discussed  above  is  not  fundamen¬ 
tally  a  relativistic  effect.  Thus  if  the  electron  energy  V 
is  much  less  than  511  keV,  so  that,  we  approximate  y^  -  1  = 
(2eV/mc^) ,  we  can  write  Eq.  (5)  as 


(VBqJct  =  [(Stt 2mV/eB2*2)1/3  -i]372 


[(21.2V1/2/B  i  )2/3  -  1]3/2 
o  o 


(6) 


where,  in  the  final  expression,  V  is  in  volts,  BQ  is  in 
gauss,  and  Z0  is  in  cm. 

In  the  experiments  to  be  described,  electron  beams  in 
the  energy  range  4-14  keV  were  employed;  a  simple  dc  low- 
current  </v  10 ?s  of  mA)  crt  electron  gun  could  then  be  used 
to  provide  the  electron  beam  with  a  diameter  of  about  1  mm 
and  energy  resolution  of  better  than  1%.  The  helical  wig¬ 
gler,  to  be  described  more  fully  below,  had  a  period  lQ  = 
3.6  cm.  Thus,  from  Eq.  (6),  one  sees  that  the  transition 
from  stable  to  unstable  orbits  would  occur  for  very  small 
wiggler  fields  indeed  if  the  axial  magnetic  field  were 
adjusted  to  be  slightly  above  5.89V1/2  gauss,  i.e.  in  the 
range  between  350  and  700  gauss.  The  axial  magnetic  field 
was  in  fact  adjusted  to  dc  values  between  about  300  to  3000 
gauss.  For  a  given  electron  energy  V  and  axial  field  B0 
the  wiggler  field  amplitude  B_u  was  varied  continuously  in 
time  by  triggering  a  spark  gap  to  discharge  a  capacitor  in 
series  with  the  wiggler  coil.  The  ensuing  RC-decay  could 
be  calibrated  to  give  B+  values  as  a  function  of  time 
during  each  discharge  puise. 

The  wiggler  coil  itself  was  a  bifilar  periodic  winding 
of  3  ram  diam  conductor  would  on  a  53  mm  diam  cylinder  with 
a  uniform  pitch  of  36  mm.  The  uniform  portion  was  666  mm 
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long,  i.e.  18.5  periods.  At  each  end  the  wiggler  diameter 
tapered  outward  to  100  mm  over  a  175  mm  length.  It  was 
found  that,  in  addition  to  provision  of  these  tapered  end 
portions,  careful  symmetrizing  of  the  conductors  at  the  end 
turns  was  essential  for  obtaining  stable  beam  transmission 
through  the  wiggler.  Furthermore,  flux  shunts  at  the  ends 
were  required  to  produce  a  smooth  uniformly  tapered  transi¬ 
tion  into  the  wiggler.  A  plot  of  one  component  of  the 
transverse  field  produced  by  this  wiggler  is  shown  in  Fig. 
2.  (The  uniform  portion  is  not  shown,  as  this  portion  is 


Figure  2.  Measured  transverse  magnetic  field 
at  the  entrance  end  of  the  wiggler. 


relatively  easy  to  produce.)  This  wiggler  produced  a  field 
of  about  20  gauss/kA,  and  fields  up  to  250  gauss  have  been 
routinely  produced. 

Several  beam  analyzers  were  constructed  to  examine  the 
properties  of  the  beam  within  the  uniform  portion  of  the 
wiggler.  For  the  data  to  be  presented  in  this  paper,  a 
movable  analyzer. was  used  consisting  of  two  parallel  plates 
spaced  by  9  mm  and  positioned  normal  to  the  axial  magnetic 
field.  The  first  plate  had  a  3  mm  hole  in  its  center 
through  which  the  beam  would  pass  either  in  the  absence  of 
any  wiggler  field,  or  for  wiggler  field  values  below  the 
critical  value.  In  this  case,  paraxial  helical  orbits  with 
diameter  less  than  3  mm  were  ascertained  to  be  produced,  so 
that  the  beam  current  was  collected  by  the  back  plate.  If 
the  orbit  were  to  involve  excursions  of  more  than  3  mm  away 
from  the  axis,  current  would  be  collected  by  the  front 
plate.  When  the  beam  was  seen  to  migrate  back  and  forth 
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between  the  two  analyzer  plates  as  the  wiggler  field  decayed 
with  time,  this  was  taken  as  direct  evidence  for  a  strongly 
non-helical  orbit.  Two  examples  of  this  migration  are 
shown  in  Fig.  3,  which  is  traced  from  oscillograms  of  the 


TIME  (msec) 


Figure  3.  Measured  currents  to  front  (f) 
and  back  (b)  plates  of  beam  analyzer. 
Arrows  indicate  abrupt  transitions  from 
non-helical  to  helical  orbits.  Lower 
example  is  for  a  lower  axial  field  value 
than  upper  example,  so  that  transition 
occurs  at  higher  value  of  wiggler  field. 


current  waveforms  to  the  analyzer  plates  as  a  function  of 
time  following  firing  of  the  wiggler  field  spark  gap.  The 
examples  are  for  two  different  axial  field  values  (lower 
for  the  bottom  example  than  for  the  top) .  One  sees  the 
beam  gyrate  wildly  back  and  forth  between  the  two  plates 
until  a  certain  time,  denoted  by  the  arrows,  when  the  wig¬ 
gler  field  has  decayed  to  a  specific  value.  The  transition 
to  beam  collection  by  the  back  plate  alone  (i.e.  paraxial 
helical  orbits)  is  seen  to  be  abrupt.  Values  of  wiggler 
field  were  noted  at  each  transition  point  observed  when 
axial  field  and  beam  energy  were  varied.  These  values  are 
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plotted  in  Fig.  4  as  a  function  of  the  independent  variable 


Figure  4.  Measured  values  of  B.u/B0  at  which 
transitions  from  stable  to  unstable  orbits 
were  observed,  for  electron  energies  between 
4-14  keV.  Solid  line  is  theoretical  pre¬ 
diction. 

(S^V^/Bq)2/3  -  1,  as  suggested  by  Eq.  (6)  for  l0  * 

3.6  cm.  The  straight  line  in  Fig.  4  is  this  same  variable 
raised  to  the  three-halves  power. 

Transitions  from  stable  to  unstable  orbits  have  been 
observed  for  wiggler  fields  as  low  as  2  gauss  (lowest 
datum  in  Fig,  4). 


Discussion 


Magnetic  wigglers  for  free  electron  laser  applications 
produce  helical  electron  orbits  in  the  absence  of  an  axial 
guide  field,  but  may  produce  strongly  non-helical  orbits 
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if  an  axial  field  is  present.  One  predicted  (8)  conse¬ 
quence  of  this  phenomenon  is  an  abrupt  jump  in  the  orbit 
from  non-helical  to  helical  once  the  magnetic  wiggler  field 
strength  falls  below  a  critical  value,  for  fixed  axial 
field  and  beam  energy.  This  behavior  has  been  observed  ex¬ 
perimentally  over  a  wide  range  of  (non-relativistic)  beam 
energies  and  axial  field  strengths.  The  data  follow  an 
approximate  three-halves  power  law  in  the  variable 

*■'  -  1,  as  suggested  by  the  theory.  The  data 

fall  systematically  about  10-20%  higher  in  this  variable 
than  is  predicted  (corresponding  to  about  a  factor-of-two 
smaller  value  of  B4-/BQ  than  is  predicted) .  An  overestimate 
in  measured  electron  beam. energies  could  explain  the  dis¬ 
crepancy  between  theory  and  experiment,  but  measurement 
accuracies  are  believed  sufficient  to.  rule  this  out.  Finite 
geometry  effects,  due  either  to  off-axis  departures  of  the 
wiggler  field  from  Eq.  (1),  or  from  the  finite  spatial 
resolution  of  the  analyzer,  could  also  contribute  to  the 
apparent  discrepancy. 

However,  the  crucial  points  for  users  of  magnetic  wig- 
glers  in-  axial  guide  magnetic  fields  are  (1) ,  the  care  re¬ 
quired  in  wiggler  construction  (especially  at  the  "first" 
turn,  and  within  a  gradual  transition  region)  in  order  to 
observe  a  paraxial  helical  orbit  at  all;  and  (2) ,  the  clear 
observation  of  an  abrupt  transition  between  stable  and  un¬ 
stable  orbits  at  (sometimes  very  low)  critical  wiggler 
fields,  much  as  had  been  predicted  by  theory. 

It  may  be  that  the  non-helical  orbits  will  be  of  util¬ 
ity,  although  it  would  be  easy  to  despair  in  attempting  to 
formulate  a  theory  for  free  electron  laser  operation  with 
such  a  complex  equilibrium  state.  These  orbits  can  possess 
large  amplitude  harmonic  overtones  (10)  which  should  radiate 
incoherent  radiation  at  wavelengths  a  few  time  shorter  than 
l0/2y^ .  It  may  even  be  possible  to  observe  coherent  ampli¬ 
fication  on  such  a  spatial  overtone  of  the  fundamental 
wiggler  period;  but  speculation  carries  risks.... 
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DEGRADATION  IN  GAIN  FOR  A  FREE  ELECTRON  LASER 
AMPLIFIER  DUE  TO  ELECTRON  MOMENTUM  SPREAD* 


A. 
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t 


Center  for  Plasma  Physics 
and  Racah  Institute  of  Physics 
Hebrew  University  of  Jerusalem,  Israel 


A  finite  spread  in  axial  momentum  for  the  electron  beam  in 
a  free  electron  laser  amplifier  is  shown  to  decrease  the 
small-signal  gain.  For  millimeter  and  sub-millimeter  wave 
amplifiers,  where  exponential  growth  dominates  the  gain,  it 
is  shown  that  the  gain  is  approximately  3  db  below  that  for 
a  cold  beam  if  the  relative  momentum  spread  (&u/u)i/2  = 
(GQ/248)  l'2(\0/L)  ,  where  G0  >>  1  is  the  gain  in  db  for  the 
cold-beam  case,  A0  is  the  magnetic  wiggler  period,  and  L  is 
the  amplifier  length.  Exact  numerical  examples  are  given 
for  representative  FEL  amplifiers  at  35  and  550  GHz. 

Key  words:  free  electron  laser,  amplifier,  electron  momen¬ 
tum  spread. 


Most  theoretical  work  concerning  amplification  of 
radiation  in  free  electron  lasers  (FELs)  deals  of  necessity 
with  idealized  models.  One  idealization  widely  employed 
involves  the  neglect  of  finite  momentum  spread  of  the  elec¬ 
tron  beam.  The  underlying  mechanism  for  small-signal  ampli¬ 
fication  involves  axial  synchronization  in  propagation  ve¬ 
locity  between  one  of  the  allowed  modes  of  radiation  sup¬ 
ported  by  the  beam,  and  the  beam  itself.  Thus  when  a  spread 
in  axial  beam  momentum  is  present,  a  mixing-in-phase  can  be 
expected  to  degrade  the  amplification  which  would  otherwise 
be  predicted  for  a  cold  beam.  Prior  workers  (1,2)  have 
taken  note  of  this  fact  and  have  provided  estimates  of  the 
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effect  of  momentum  spread.  This  paper  presents  an  exact 
analytical  model  to  account  for  finite  momentum  spread  for 
a  particular  distribution  function.  When  exponential 
growth  dominates  the  gain,  a  simple  approximate  formula  is 
derived  to  estimate  the  loss  in  gain  due  to  the  momentum 
spread.  Exact  numerical  examples  are  also  given  for  repre¬ 
sentative  FEL  amplifiers  at  35  and  550  GHz. 

The  basic  FEL  model  adopted  here  is  identical  to  that 
treated  by  Bernstein  and  Hirshfield  (B-H)  (3).  That  work 
gave  an  exact  small-signal  solution  of  the  Vlasov-Maxwell 
equations  for  the  steady-state  evolution  of  the  co-propa- 
gating  disturbance  which  grows  in  space  on  a  relativistic 
electron  beam  passing  along  the  axis  of  a  helical  magnetic 
wiggler.  The  B-H  theory  was  derived  for  a  beam  of  arbitra¬ 
ry  momentum  distribution  in  a  wiggler  of  arbitrary  strength, 
but  the  solutions  presented  were  for  the  case  of  a  cold 
beam,  viz., 

f0(a,B,u)  =  No6(a)6(6)6(u-U)  .  (1) 

Here  a  and  0  are  the  two  transverse  components  of  canonical 
angular  momentum  Ux  -  eAx/mc2  and  Uy  -  eAy/mc2,  Ax  and  Ay 
are  the  components  of  the  wiggler fs  vector  potential,  Ux 
and  Uy  are  the  transverse  components  of  translational  momen¬ 
tum,  and  U  *  (y2  -  l)1/2  is  the  total  momentum  as  related 
to  the  relativistic  energy  factor  .  (All  momenta  are 
normalized  to  me.)  Eq.  (1)  thus  describes  a  beam  which, 
prior  to  entering  the  wiggler,  contains  electrons  possess¬ 
ing  both  zero  transverse  momentum  and  unique  axial  momen¬ 
tum  U. 


As  mentioned  above,  an  important  source  of  degraded 
amplification  is  the  finite  spread  of  axial  momentum  on 
the  electron  beam.  In  the  work  reported  here,  we  choose 
the  simplest  distribution  capable  of  describing  such  a 
spread,  viz., 


fQ(a,6,u) 


NQ6(a)6(B) 


^H(u-U^)  -  H(u-U2) 
AU 


(2) 


where  H(x)  =  1  for  x  >  0,  H(x)  *  0  for  x  <  0,  and  AU  = 

U2  -  Ui  >  0.  This  distribution  can  of  course  not  be  real¬ 
ized  in  nature  [in  the  same  sense  that  the  distribution 
given  by  Eq.  (1)  cannot].  It  may,  however,  not  be  a  bad 
approximation  for  certain  accelerators  (except  for  the 
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sharp  edges) ;  but  its  utility  here  is  that  it  enables  an 
analytic  form  to  be  derived  for  the  governing  dispersion 
relation. 


The  goal  of  the  present  work  is  identical  to  that  in 
B-H,  namely  to  calculate  the  power  gain  G  (in  db)  for  a 
single  pass  of  electromagnetic  radiation  along  a  FEL  ampli¬ 
fier  of  length  L. 

G/10  * 

10W1U  =  a2(L)a2(L)  -  1  .  (3) 


Here  a2(L)  is  the  dimensionless  wave  electric  field  at  the 
amplifier  output,  normalized  to  unity  at  the  input.  The 
subscript  ,,2M  labels  one  of  the  three  polarizations  per¬ 
mitted,  namely  that  which  twists  in  space  a  quarter-period 
behind  the  wiggler's  vector  potential.  [Eqs.  (35)  and  (37) 
in  B-H  give  the  other  two  polarizations.) 


The  wave  amplitude  a2(L)  Is  a  superposition  of  several 
co-propagating  normal  modes,  each  with  its  wavenumber  k j , 
viz.  , 


a2(L) 


R’OCj) 


exp(ikjL) 


(4) 


The  relative  mode  amplitudes  B(kj)/R* (kj)  are  prescribed 
once  boundary  conditions  are  set.  R(kj)  *  0  is  the  disper¬ 
sion  relation  for  the  system  which  determines  the  kj (u) , 
assuming  R_1(k)  to  have  simple  poles.  For  the  cold  beam 
case  R(kj)  is  a  sixth-order  polynomial. 


R(x) 


[(x-p)2  -  62(1+£2)][(x+xo)2  -  b2][(x-xo)2 

,  ,.2.2,  2  .2.  ,  2.2  ,2. 

+  £  6  (x  -b  ) (x  +x0~b  )  , 


where  x  *  kc/w,  xQ  *  k0c/w;  6  •  (top/w)  (U/yU|)  ^2,  b  = 
(1-U262)1/2,  v  -  Y/Uz,  Uz  =  (U2-£2)l/2,  and  £  *  -  eB0/mc2k0. 
The  wiggler  field  strength  and  wavenumber  are  B0  and  k0. 

This  equation  has  been  obtained  as  well  by  Sprangle  (1), 
and  related  forms  have  been  derived  and  discussed  by 
Kroll  and  McMullin  (2)  and  by  Kwan,  Dawson,  and  Lin  (4). 

When  6  «  Xo  <<  1»  a  reduced  form  of  Eq.  (5)  is  a  good 
approximation,  namely 
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R(x)  -  [(x-y)2  -  62(1+C2)][x  -  (b+x  )]  +  4  5262x  .  (6) 

o  l  o 

For  kQ/k  £  (1+C2)/2y2  the  maximum  growth  occurs  near  b+x0 
-  U  =  (5252x0/2) .  To  requisite  accuracy  the  roots  are 

x^  *  y  +  (6^^^xo/2)^^exp(-7ri/3) 


/  r2  2  /ov  1/3 
x,  s  p  -  (5  5  x  /2) 
o  o 

These  roots  are  of  use  in  scaling  estimates  when  exponen¬ 
tial  gain  is  dominant.  Exact  numerical  evaluations  given 
in  B-H  show,  however,  that  Eq.  (7)  cannot  be  used  to  deter¬ 
mine  the  entire  gain  spectrum. 

When  Eq.  (2)  is  employed  as  the  distribution  function 
all  the  momentum-space  integrals  in  the  Vlasov  formulation 
can  be  expressed  analytically.  We  then  find 


R(x)  =  [ (x-y1) (x-y2)  -  6 ' 2(l+£2) ] [ (x+xq) 2  -  b * 2 ] 

r/  \2  vt2,  ^  r2.,2,  2  ,2W  2.  2  ,  ,2N 

x  [(x-xq)  -  b  ]  +  £  6*  (x  -d  )  (x  -rxQ-b  ) 

where 


2 


Au  =■  U2  -  Vi  <  0,  Y?,2  "  1  +  U2>2,  Ul1>2  -  U2>2  -  C2, 

hi  ,2  *  Yi  t2^Zl  ,2  • 


(8) 
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and 
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When  AU/U  «  1,  Ay  *  -  UAU(1  +C2)yU|,  6'  =  6,  and 
b*  s  d  -  b.  Thus  the  only  effect  of  finite  momentum  spread 
in  this  limit  is  in  the  factor  (x  -  Ui)(x  -  uO  *  (x  -  u)2 
-  (Ay/2)2  in  the  first  bracket  in  Eq.  (8),  where 
y  =  (y i  +  U2 ) /2 .  The  close  similarity  between  Eqs.  (8)  and 
(5) ,  and  the  simplicity  of  the  former,  make  determination 
of  the  roots  kj  a  routine  matter.  This  simplicity  is  not 
enjoyed  when  the  momentum  spread  is  described  by  functions 
U)  with  non-zero  values  of  3f0/3U  in  a  finite 
interval,  because  of  wave-particle  resonance  effects. 

As  for  the  cold-beam  case,  where  6*  <<  Xq  <<  1,  Eq. 

(8)  may  be  reduced  to  the  approximate  form 


R(x)  =  [(x-U)2  -  (Ay/2) 2]  [x  -  (b'+xQ)]  +  =  0.  (9) 


If  (Ay/2)2  «  3(§25’2x0/2)1/3,  the  roots  of  Eq.  (9)  near 
b'  +  x0  -  y  =  (C26’2xo/ 2)1/2  are  approximately 

-  7  1/7 

=  y  +  (5  6  *  *~Xq/2)  exp (i*rr/3) 

+  |(Ay/2)2U2«,2xo/2)"1/3e:;?(-iT/3) 

(10) 

/r2.,2  ...1/3  1,.  ,..2.. 2. ,2  ,-.-1/3 

x,  =  y  -  (C  <5  x  /2)  -  -r(Ay/2)  (5  6'  x  /2) 

J  o  J  o 


Thus  the  spatial  growth  constant  Inuq  is  seen  to  decrease 
on  account  of  momentum  spread  as 


(ID 


For  pure  exponential  gain,  i.e.  excluding  the  15.6  db  input 
coupling  loss  (see  B-H) ,  one  has 


G  =  54.58(L/A)Imx1  db  (12) 

where  X  is  the  radiation  wavelength.  From  Eq.  (11)  we  can 
write  G  *  G0  -  Gj,  where  G0  is  the  gain  with  no  momentum 
spread,  and  G^  is  the  small  decrease  due  to  the  spread 
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G  =  54.58a/A)^  U2S2x  /2)1/3  db  .  (13) 

o  z  o 

For  £  =  0.47,  X  -  4.9,  x0  =  2.73  *  10~2 ,  62  =  3.80  x  1CT6 , 
and  L/A  «  367  (corresponding  to  a  representative  FEL  ampli¬ 
fier  to  be  discussed  below),  Eq.  (13)  gives  G0  *  39.1  db . 
[If  one  subtracts  the  15.6  db  input  coupling  loss,  the 
actual  gain  would  be  23.5  db  (at  a  wavelength  of  560  ym)  .  ] 
Now 


54.58  L  n  2,1  _2 *2  .-1/3  .. 

-  Y  (t  5  5  *  )  db  •  (14) 

8/j  A  2  o 


Substituting  from  Eq.  (13)  gives  the  value  of  Ay  which 
would  bring  about  a  gain  loss  Gj 

(Ay)2  =  5.37  x  10_3G  G. (A/L)2  .  (15) 

o  1 

For  the  example  cited  above  with  L/A  *  367  we  find  Ay  = 

2,16  x  10  3  for  G0  =  39.1  db  and  =  3  db ,  i.e.  for  a 
factor-of-two  decrease  in  power  amplification.  This  corre¬ 
sponds  to  a  relative  momentum  spread  AU/U  * 

| Ay  1 [yU|/U2(l  +  C2)]  of  0.041. 

Equation  (10)  also  suggests  that  the  frequency  at 
which  gain  has  its  peak  value  will  decrease  as  momentum 
spread  increases. 


Exact  numerical  evaluations  for  small-signal  gain  G 
have  been  carried  out  using  the  full  dispersion  relation 
[Eq.  (8)],  and  with  amplitudes  [see  Eq.  (4)]  appropriate 
to  a  perfectly  matched  amplifier  output.  One  example  is 
for  a  mm-wave  amplifier  employing  an  electron  beam  typical 
of  that  produced  by  a  small  Febetron  accelerator,  with  y  * 
1.78,  J  =  100  A/cm2,  AQ  -  3.6  cm,  £  =  0.2,  and  L  =  36  cm. 
Gain  curves  are  shown  in  Fig.  1  for  zero  momentum  spread, 
and  for  finite  momentum  spreads  between  5  and  20%.  Gain  is 
seen  to  fall  by  one-half  for  AU/U  =  0.15,  and  the  frequency 
for  peak  gain  drops  by  about  6%.  A  second  example  is  for  a 
sub-mm  wave  amplifier  employing  a  beam  typical  of  the  VEBA 
accelerator  at  Naval  Research  Laboratory,  with  y  *  4.9, 

J  =  6  kA/cm2,  A0  =  2.0  cm,  £  =  0.47,  and  L  =  20  cm.  For 
this  case  the  computed  gain  curves  are  shown  in  Fig.  2, 
again  for  zero  momentum  spread  and  for  spreads  between 
5  and  20%. 
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Figure  1.  Gain  curves  for  a  FEL  using  a  400  kV 
electron  beam,  for  electron  momentum  spread 
between  0  and  20%. 


Frequency  (GHz) 

Figure  2.  Gain  curves  for  a  FEL  using  a  2.0  MV 
electron  beam,  for  electron  momentum  spread 
between  0  and  20%. 
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Comparisons  between  the  exact  results  (Fig.  2)  and  the  * 
approximate  predictions  [Eqs.  (12-15)]  are  instructive. 

The  peak  gain  for  the  cold  beam  is  17.8  db  (i.e.  60*)  com¬ 
pared  with  the  approximate  value  of  23.5  db.  The  gain 
drops  by  half  to  14.8  db  (i.e.  30x)  for  AU/U  somewhat 
greater  than  10%;  our  approximate  result  is  4.1%.  These 
comparisons  for  the  example  presented  in  Fig.  1  are  not 
meaningful  since  the  peak  gain  G0  is  less  than  7.8  db  (6x) . 

Finally,  we  point  out  the  scaling  laws  suggested  by 
Eqs.  (12-15),  valid  for  high  gain  devices  where  exponential 
growth  dominates.  For  negligible  momentum  spread. 


G  „  J^V'V2'3*-1'3 


db 


or  equivalently 


G 

o 


•v  j1/3(l/»o){2/V/3 


db 


(16) 


For  the  gain  decrease  Gx  «  G0  due  to  finite  momentum 
spread  we  have 


G1Gq  -v  (AU/U)2(L/Aq)2  (db)2  .  (17) 

Eq*  (17)  indicates  that  high  gain  short  amplifiers  are  less 
susceptible  to  gain  degradation  due  to  momentum  spread, 
than  are  low  gain  long  amplifiers.  This  scaling  is  inde¬ 
pendent  of  A  and  y  provided  G0  is  high.  For  Gj  *  3  db,  the 
numerical  value  for  Eq.  (17)  gives  (AU/U)i/2  =  (Gq/248)1/2 
(A0/L),  where  (A?J/U)j/2  relative  momentum  spread  for 

a  factor-of-two  decrease  in  gain.  Gain  degradation  for 
long-wiggler  FELs  operating  in  the  collective  regime  can 
be  expected  to  be  serious  unless  AU/U  <<  1. 

It  should  be  added  as  a  caveat  however  that  momentum 
spread  may  not  always  degrade  gain  in  a  FEL.  The  geomet¬ 
rical  optics  theory  for  a  FEL  amplifier  (5)  shows  that 
gain  may  arise  from  a  wave-particle  resonance,  provided 
fQ(a,6>u)  is  not  symmetric  in  u  about  its  maximum,  and 
provided  df0/du  has  the  requisite  sign  at  the  wave’s  phase 
velocity.  It  is  expected  that  this  mechanism  would  compete 
with  that  discussed  in  the  present  paper,  and  could  in  fact 
allow  substantial  gain  in  the  presence  of  tailored  momentum 
spread. 
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